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Abstract 

We study the four-point correlation function of stress-tensor supermultiplets in Af = 4 
SYM using the method of Lagrangian insertions. We argue that, as a corollary of M = 4 su- 
perconformal symmetry, the resulting all-loop integrand possesses an unexpected complete 
symmetry under the exchange of the four external and all the internal (integration) points. 
This alone allows us to predict the integrand of the three-loop correlation function up to 
four undetermined constants. Further, exploiting the conjectured amplitude/correlation 
function duality, we are able to fully determine the three-loop integrand in the planar 
limit. We perform an independent check of this result by verifying that it is consistent 
with the operator product expansion, in particular that it correctly reproduces the three- 
loop anomalous dimension of the Konishi operator. As a byproduct of our study, we also 
obtain the three-point function of two half-BPS operators and one Konishi operator at 
three-loop level. We use the same technique to work out a compact form for the four-loop 
four-point integrand and discuss the generalisation to higher loops. 
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1 Introduction 



In the first few years after the discovery of the AdS / CFT correspondence [T] considerable 
effort was put into studying gauge-invariant operators and their correlation functions in 
perturbative M = 4 SYM, see for example [2T4TT]. 

In parallel a great deal of progress has been made in exploring seemingly completely 
different objects in M = 4 SYM, namely scattering amplitudes. One of the important 
results in this arena has been the discovery that planar MHV amplitudes can be described 
by null polygonal Wilson loops [T2HT4] . This has led to much progress in understanding the 
amplitudes themselves (since so far the Wilson loop integrals have been easier to compute 
than the corresponding amplitude integrals [T5TT2T| ). culminating in the discovery of a 
new symmetry, dual conformal symmetry [IB] and dual superconformal symmetry [2"2"H2"4] . 
This in turn has been instrumental in recent progress in studying scattering amplitudes 
at the level of the integrand beyond the MHV case j25j[26] and in the construction of a 
supersymmetric Wilson loop [271128]. which should reproduce the super-amplitude in planar 
Af = 4 SYM provided the correct regularisation is understood [2H] . 

Excitingly, a year ago it was discovered that correlation functions of gauge-invariant 
operators, defined as superconformal primaries of M = 4 supermultiplets, describe both 
Wilson loops [30] and MHV amplitudes [31]. These new dualities require taking the light- 
like (or "on-shell") limit in which the operators are located at the vertices of a null poly- 
gon. The correlation function/Wilson loop duality can be viewed as a relation between the 
(properly regularised) integrals in terms of which both objects are expressed. However, 
the correlation function/amplitude duality is most easily understood at the level of the 
integrand [32]. The natural definition of the integrand for loop- level correlation functions 
via the Lagrangian insertion procedure remarkably reproduces the amplitude integrand 
derived in momentum twistor space [26]. Both the MHV amplitude and the correlation 
function of superconformal primaries have natural superspace generalisations, and indeed 
(if we switch off the anti-chiral half of the superspace coordinates), the full supercorrela- 
tion functions of the supermultiplets continue to reproduce all non-MHV superamplitude 



We expect that this duality will lead to new insights both for correlation functions 
and amplitudes, and it is the purpose of this paper to give an example. Namely, we shall 
exploit the duality, combined with a new symmetry to be described below, to derive the 
correlation function of four stress-tensor supermultiplet operators at three and four loops 
and to set up a procedure for going to higher loops. 

The stress-tensor multiplet operator T is the simplest example of the so-called half-BPS 
multiplets. Its lowest-weight state (superconformal primary) is a scalar operator O 2 o' °f 
scaling dimension two in the 20' of SU(4). In the context of the AdS/CFT correspondence 

1 We should note that it is expected that this duality holds for correlation functions of a large class 
of operators, but so far has only been checked for the stress-tensor multiplet and the Konishi 

multiplet |30}135| . An explanation for this duality was proposed in momentum twistor space |35) . However 
this explanation is at a rather formal level due to the regularisation issue alluded to earlier, so at the 
moment we view the duality as conjectural. 
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the half-BPS operators are regarded as the duals of the AdS 5 x S 5 supergravity states. 
They have the important property that their two- and three-point functions are protected 
from quantum corrections. The first non-trivial example of a correlation function of such 
multiplets is provided by the four-point case. In the past the four-point functions of the 
operators C20' have been extensively studied, both at strong coupling (in AdS supergravity) 
[3o] and at weak coupling [2"HTTj. The perturbative calculations used the standard Feynman 
diagram technique adapted to Af = 1 superspace [I] or to N — 2 harmonic superspace [3]. 
Beyond two loops this technique becomes exceedingly difficult, so no explicit three- loop 
results for four-point functions are available. Our aim in this paper is not to carry out a 
Feynman diagram calculation, but to predict the form of the correlation function in terms 
of a set of scalar Feynman integrals. 

There are two main ingredients which allow us to do this. Firstly, we discover a hid- 
den permutation symmetry of the integrand of the four-point correlation function, which 
exchanges the integration points and the external points. This symmetry can be seen 
in the framework of the Lagrangian insertion procedure. The I— loop correction to the 
four-point correlation function of the scalar operators O20 1 is given by their (4 + €}— point 
Born-level correlation function with i Lagrangian insertions. The M = 4 SYM Lagrangian 
is a member of the same stress-tensor super mult iplet T as the operator O20' itself. Con- 
sequently, what determines the integrand of the four-point I— loop correlation function 
is the Born-level correlation function of (4 + £) operators T . It can be described by a 
nilpotent superconformally covariant (4 + £)— point polynomial in M = 4 analytic super- 
space [SJG2]. Our new observation is that this nilpotent object has full permutation S^+e 
symmetry. This, together with the crossing symmetry of the correlation function of the 
super-operators T, lead to the aforementioned new symmetry of the four-point integrand. 

This hidden symmetry, together with reasonable assumptions about the pole structure 
of the integrand following from the OPE, already yield a very constrained form for the 
four-point correlation function integrand at any loop level. At three loops, it determines 
the result, for an arbitrary gauge group (e.g., SU(N C ) for any value of N c ), up to four 
unknown coefficients. 

The second main ingredient in our construction is the correlation function/amplitude 
duality. According to it, the integrand of the planar four-point correlation function of 
protected operators O 2 o', in the light-like limit where the four points become sequentially 
light-like separated, is equal to the square of the integrand of the planar four-gluon scat- 
tering amplitude. We use this duality to compare the most general form of the four-point 
correlation function described above with the known expression for the three-loop four- 
gluon amplitude [38]. We find that, for a certain choice of the four coefficients, we get 
precise agreement, thus at the same time determining the unique form of the planar corre- 
lation function at three loops and providing further non-trivial evidence for the duality. It 
is interesting to notice that, besides the known three- loop integrals of "ladder" and "tennis 
court" type, we find only two new three-loop planar integral topologies. 

Our procedure is not limited to three loops. It is rather straightforward to use the 
results on the four-gluon amplitude at four [39] and even five [40] loops and predict the 
form of the integrand of the correlation function. In this paper we treat the four-loop case 
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and obtain a compact expression for the correlation function integrand there. 

Although in this paper our focus is on the correlation functions, we wish to emphasise 
that the techniques have much to say about the amplitude integrands as well. In particular 
the new permutation symmetry, which relates external points and integration points, is 
also present for the amplitudes (albeit in a broken form). This together with the fact that 
the amplitude/correlation function duality relates products of lower-loop amplitudes to 
the correlation function means that a large portion of the amplitude can be found from 
lower-loop amplitudes. For example, the entire three-loop amplitude can be determined 
from the one- and two-loop amplitudes and a large portion of the four-loop amplitude as 
well. We comment further on this point in the conclusions. 

The paper is organised as follows. In Section [2] we define the correlation functions of the 
stress-tensor multiplet and give a short review of the Lagrangian insertion technique. In 
Section [3] we derive the three-loop correlation function and check the agreement with the 
four-gluon scattering amplitude in the light-like limit. In Section H] we test the correlation 
function by an OPE analysis and we show that it correctly reproduces the known value 
for the three-loop Konishi anomalous dimension j?T]. This section also contains our pre- 
diction for the three-point correlation function of two protected scalar operators and the 
Konishi operator at three loops. In Section we derive a similar expression at four loops 
albeit without the benefit of the analogous OPE check. In the conclusions we comment 
on the implications of these techniques for amplitudes and also discuss the conditions for 
planarity. In Appendix [B] we prove the hidden permutation symmetry which is the key to 
the simple form we find for the three-loop correlation function. It is due to the existence of 
an S^— invariant, conformally covariant, nilpotent polynomial which we give a simple ex- 
pression for here. Appendix [C] contains some technical details on the numerical evaluation 
of the integrals needed for the OPE test. 

2 Correlation functions of the stress-tensor multiplet 
in AT = 4 SYM 

In this section, we define the correlation functions of the stress-tensor multiplet in the 
M = 4 SYM theory and discuss their general form. We begin by presenting the notation 
(which is inherited from [SSIEI]) and summarise some basic results concerning correlation 
functions in M = 4 SYM, following from M = 4 superconformal symmetry. 

2.1 The J\f = 4 stress-tensor multiplet in analytic superspace 

The field content of the M = 4 super- Yang-Mills theory comprises six real scalars $ 7 (with 
I = 1, ... ,6 being an SO(6) index), four Weyl fermions (and their complex conjugates) 
and the gauge field, all in the adjoint representation of, e.g., the SU(N C ) gauge group. In 
this paper, we study the correlation functions of the special class of the so-called half-BPS 
gauge-invariant local operators. One of the important properties of such operators is that 
their scaling dimension is protected from perturbative corrections. 
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The simplest example of a half-BPS operator is the bilinear gauge-invariant operator 
made of the six scalars, 

Oi J , = tr - ^5 IJ ti ($ A '$ K ) . (2.1) 

The operator O u belongs to the representation 20' of the R symmetry group 5*0(6) ~ 
5{7(4). It is the lowest-weight state of the so-called M = 4 stress-tensor supermultiplet 
containing, among others, the stress-tensor (hence the name) and the Lagrangian of the 
theory. The latter appear as the coefficients in the expansion of half-BPS superfield op- 
erator T(x, 9 , 9a) in powers of the odd variables 9 A and 9 \ (with A = 1, ... ,4 being an 
SU(A) index). 

The half-BPS superfield T(x, 9 A , 9 a) satisfies constraints implying that it depends on 
half of the odd variables, both chiral and anti-chiral. The appropriate formalism which 
makes the M = 4 half-BPS property manifest is that of M = 4 analytic superspace [37] or, 
equivalently, M = 4 harmonic superspace [12] . On the other hand, the lack of an off-shell 
formulation of M = 4 SYM makes M = 4 analytic superspace inadequate for performing 
Feynman graph calculations of the correlation functions. These are most conveniently done 
in M = 2 harmonic superspace [13], after which the results can be lifted to M = 4. 

Harmonic/analytic superspace makes use of auxiliary bosonic ("harmonic") y— variables 
in order to covariantly break the R symmetry group 577(4) down to SU (2) x 577 (2)' x £7(1). 
Specifically, the chiral Grassmann coordinate 9 A can be decomposed into its halves by 
splitting the 577(4) index A = (a, a') (with a, a' = 1, 2) and substituting 

0a^GCO, with p£ = 0- + #i£, (2.2) 

and similarly for the anti-chiral 9°a — > (p"/ ,9%)- Here the harmonic variables y^, carry 
577(2) x SU (2)' indices and parametrise the four- dimensional coset 577(4) / (SU(2) x SU (2)' x 
i7(l)) (or its complexification in the approach of [3Z])- The projected odd variables carry 
also a {7(1) C 5{7(4) weight (+1) for p and (—1) for p. Then, the stress-tensor supermul- 
tiplet is defined as a function on analytic superspace 

T=T(x aa jP a aJ pZ„y a al ). (2.3) 

The half-BPS nature of T manifests itself in the independence of this function on 0%' and 9%. 
This is the meaning of the so-called "BPS shortening" of the stress-tensor supermultiplet. 

The supermultiplet ( 12. 3p carries {7(1) weight (+4) and its scaling dimension is protected 
from quantum corrections. Its lowest-weight component (superconformal primary) is given 
by the bilinear scalar operator ( 12. ip 

0(x,y) = YiYjO&(x)=Y I Yjtr(<f> I $ J ) , (2.4) 

2 In this half-BPS supcrficlds radically differ from chiral supcrficlds F(x,0 A ) which depend only on the 
chiral half of the odd variables. 
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where the 50(6) indices have been projected with the 5*0(6) harmonic variable Yj = 
Yj(y), a (complex) null vector, Y 2 = YjYj = 0. It carries U(l) weight (+2) and can 
thus be expressed as a polynomial of degree two in the SU(4) harmonic variables y^, (see 
Appendix |A]) . The main subject of this paper is the four-point correlation function of the 
operators (12 .4p . 

The complete stress-tensor supermultiplet (I2.3P is obtained by acting on the lowest 
weight (12 .4p with half of the M = 4 supersymmetry generators realised in the analytic 
superspace: 

T(x, p, p, y) = exp [pi Q a a + p% Q«Q 0(x, y) . (2.5) 

By construction, this multiplet is annihilated by the other half of the supercharges, hence 
the name "half-BPS". In what follows we will drop the dependence on the anti-chiral 
variables (for the explanation see Appendix [B]) and will restrict T to its purely chiral 
sector, T(x, p, 0, y). The expansion of this object in the four odd variables p a a is very short 
and its top component is the (on-shell) chiral Lagrangian of M = 4 SYM: 

T{x, p, 0, y) = 0(x, y) + ... + (p) 4 C M=i (x) . (2.6) 

The explicit form of the complete expansion can be found in Ref. [53] • Notice that Cj^ = ^(x) 
is an SU(4) singlet (hence y— independent), since the entire U(l) weight (+4) of T is carried 
by the odd factor (p) 4 = Y[ a a Pa- This allows us to write the on-shell action of the M = 4 
theory SV=4 = / d 4 x C^ =i {x) as a Grassmann integral over the chiral half of analytic 
superspace [6] 

SV=4 = J Vx J d A p T(x, p, 0, y) . (2.7) 

We recall that the action SV=4 is y— independent and it is invariant under the full M = 4 
supersymmetry. 



2.2 Correlation functions of the M = 4 stress-tensor multiplet 

Let us consider the n— point correlation function of the analytic supermultiplets T re- 
stricted to its chiral sector, see f )2.6p . By construction, it depends on n copies of the (chiral) 
super-coordinates (xi,pi,yi). Viewed as a function of the 't Hooft coupling a = g 2 N c / '(4-7T 2 ) 
and Grassmann variables pi, the correlation function admits a double series expansion 

n— 4 oo 

G n = (T(l) . . . 7»> = E ° m G nU^ ■ ■ ■ ' n ) > ( 2 - 8 ) 

k=0 1=0 

where T(i) = T(xj, pi, 0, yi). As was explained in Refs. [33|l34"]. each G^) k is a homogeneous 
polynomial in the Grassmann variables pi (with i = 1, . . . , n) of degree 4k. It is accompa- 
nied by an additional power a k of the coupling, which takes into account the offset of the 
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perturbative level of the Born-type correlation functions with Lagrangian insertions. For 
k — the functions G^) define the £— loop corrections to the correlation function of the 
lowest components (12. 4p . The maximal Grassmann degree allowed by M = 4 superconfor- 
mal symmetry on the right-hand side of (12. 8p is not 4n, but 4(n — 4) (see Appendix IB1 for 
details). 

Let us consider the correlation function (I2.8P at four points, for n — 4. In this case, the 
sum over k reduces to a single term with k — 0, which is p— independent. In other words, 
the four-point correlation function of stress-tensor multiplets does not have a purely chiral 
nilpotent sector. This allows us to replace the superfield ( 12. 6 p by its lowest component 
( 12. 4p . so that (12. 8p becomes simply 

oo 

G 4 ={0(x 1 ,y 1 )...0(x i ,y i )) = G? (1,2,3,4) . (2.9) 

Here G? = G® are the four-point loop corrections, which will be the subject of this paper. 

We start with examining ( 12. 9p at tree level, for I = 0. At this order in the coupling, 
the correlation function ( 12. 9 p reduces to a product of free scalar propagators 

Y^Uy^^^)) = {Y{y ^ ys)) = ^ , (2.10) 

where yfj is defined in Appendix |A] and Xij = Xi — xj. In this way we obtain the tree- level 
expression for the correlation function (12. 9p 

(iV c 2 - l) 2 / 

Vl2 2/34 , Via ^24 , 2/41 2/23 \ fn ii\ 

G 4 (1,2,3,4) - -——I—— + —— + — — (2.11) 

<±y±H ) \"''12 X 34 x 13 A 24 x 41 x 23/ 

+ n 2 c - 1 / ^12 ^23 ^|4^|i + 2/12 yfi yfi^is + vhyhyhvli 

( /1_'7T~2^4 \ ^y«2 ™ 2 /-v«2 ™ 2 /y»2 , y» 2 /-v»2 ™ 2 /y»2 /y»2 , y» 2 ,y»2 

^<±/l ; \ x 12 u '23 x 34 u '41 x 24 x 34 x 13 x 13 ■ // 23 x 24 x 41 

where the first and second lines describe the disconnected and connected contributions, 
respectively. 

For computing the loop corrections G4 we employ the method of Lagrangian insertions. 
It relies on the observation that the derivatives of the four-point correlation function (12. 9 p 
with respect to the coupling constant can be expressed in terms of the five-point correlation 
function involving an additional insertions of the M = 4 SYM action, 

a^G 4 = / d 4 x 5 (0(xi, yi ) . . . 0(x 4 , y A )C M=A {x b )) . (2.12) 

The relation is very useful because it allows us to compute the one-loop correction to 
G\ from the Born-level expression for the five-point correlation function. In the same 
way, successively differentiating with respect to the coupling, we can express the I— loop 
correction to G4 in terms of the (integrated) Born-level correlation function with I insertions 
of the M = 4 SYM Lagrangian. Remarkably, since the Lagrangian is itself a member of 
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the stress-tensor multiplet, the latter correlation function can be identified as a particular 
component in the expansion of the (4 + £) —point correlation function ( 12. 8 p 

G^| P1= ... =P4=0 = (0(x l7 Vl )... 0(x 4 , y 4 )C(x 5 ) . . . £(x 4+ ,))(°)(p 5 ) 4 . . . (p 4+ ,) 4 , (2.13) 

where the superscript '(0)' indicates the Born- level approximation. 

In summary, the £— loop correction to the four-point correlation function (I2.9P can be 
expressed in terms of the Born- level (4+£) — point correlation functions G^ e . e at Grassmann 
level 4£, 

Gf (1, 2, 3, 4) = J d 4 x 5 . . . d A x 4+e (± J d A p 5 . . . d A p i+e Gfl/1, . . . , 4 + £)j . (2.14) 

Here the integration over the Grassmann variables effectively picks out only the component 
of G^. f displayed in ( I2.13p . The expression in the parentheses on the right-hand side of 
( I2.14p defines the integrand of the loop corrections to the four-point correlation function. 

The purpose of our study will be to determine the integrands G^. e . e without any Feyn- 
man graph calculations. Our approach is based on, firstly, fully exploiting the symmetries 
of the integrand. In particular, we reveal a new permutation symmetry involving all the 
(4 + £) points, which is a powerful restriction on the possible form of the integrand. Sec- 
ondly, using the recently discovered correlation functions/amplitudes duality [SUES], we 
compare our prediction for G^_ e . e with the known results on the £— loop four-gluon scat- 
tering amplitude. We demonstrate that this information is sufficient to unambiguously fix 
the form of the integrand in (12.14p . In this paper we examine in detail the case £ = 3 and 
give some preliminary results on £ = 4; the generalisation to arbitrary values of £ will be 
considered elsewhere. 

To begin with, let us recall the known results for the one- and two-loop integrands [2HU 
(i.e. the five- and six-point Born-level correlation functions of Grassmann degree 4 and 8, 
respectively) written in a particularly suggestive way: 

Gg(l,2,3,4,5) 
(3^(1,2,3,4,5,6) 

Here in the second relation the sum runs over all permutations (a±, . . . , eg) of the indices 
(1, 2, . . . , 6) and the additional factor 1/48 ensures that we do not count the same term 
more than once. In relations (I2.15P and (I2.16p . the dependence on the Grassmann variables 
is hidden in X n (with n = 5,6), which is a nilpotent superconformally covariant function 
of all the superspace points (xi, pi, yi). It is described in detail in Appendix [Bj but for the 
moment we just formulate its three main properties needed for our purposes. 

Firstly, X n is a homogeneous polynomial in the Grassmann variables pi, . . . , p n of degree 
4(n— 4). We will only need the explicit expression for the coefficient in front of (P5) 4 . . . (p n ) 4 



2(iV c 2 -l) 

(_ 4vr 2)5 

2(iV c 2 -l) 

(-47T 2 ) 6 



xL x 



XltX 



1 

ni<i<j<5 X % 

\ A /v>2 ,y»2 /v>2 

48 Z^o-eSs o-i (72 ^ 0-3 0-4 o- 5 o- 6 



III 



2 

<i<j<6 X ij 



(2.15) 
(2.16) 
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since this is the only term which contributes to the superspace integral in (I2.14p . It is 
obtained by setting p x — . . . — p± = 0: 



I n \p 1 =..- Pi =o — ( x r2 x r3 x i4 x 23 x 24 x 34) x ^(1) 2, 3, 4) x (ps) 4 . . . (pn) 4 • (2-17) 

Here _R(1,2,3,4) is a rational function of the bosonic (space-time and harmonic) coordi- 
nates at the four external points 1, 2, 3, 4: 

d(i 9ui- y^yhyuVu > 2 2 _ 2 2 _ 2 2 \ 

-tX^X, ^, O, <±y — 2 2 2 2 V" t '13 J '24 x 12" t '34 x 14 x 23/ 
x 12 x 23 x 34 x 14 

^12^13^24^34 / 2 2 _ 2 2 _ 2 2 n 
+ ™2 ™2 „2 „2 l X 14 x 23 x 12 x 34 X 13 X 24J 
x 12 x 13 x 24 x 34 

2/l3?/l42/232/24 / 2 2 22 2 2 \ 

_l_ I cy» iT" T" 'v» 'T* 1" I 

+ ~2 ™2 „2 ™2 l X 12 x 34 x 14 x 23 X 13 X 24J 
X 13 X 14 X 23 X24 

4 4 4 4 4 4 

, 2/122/34 , 2/132/24 , 2/142/23 /n 1C \ 

H 2 2""^ 2 ST" 1 2 2~ • \ Z - L ° 

x 12 x 34 x 13 x 24 x 14 x 23 



It is easy to check that R(l, 2, 3, 4) as well as the first factor on the right-hand side of (I2.17P 
have £4 symmetry under the permutation of the four points. Similarly, the product of odd 
variables in (I2.17P has an obvious S* n _4 symmetry under the exchange of the remaining 
(n — 4) internal (integration) points, so that the total symmetry of (I2.17P is S4 x S , n _ 4 . 

The second crucial feature of X n is that the complete expression for it as a function 
of the n superspace points (xi,pi,yi) is invariant under exchange of any two points, i.e. 
it is completely S n permutation symmetric. For the particular component of X n given in 
(I2.17P this symmetry is broken down to 5*4 x SV1-4, but the S n symmetry is recovered after 
restoring the dependence on p 1; . . . ,p 4 (see Appendix IB]) . 

The third feature of X n is that it has (the chiral half of) M = 4 superconformal symme- 
try. Namely, it is invariant under chiral super-Poincare symmetry and anti-chiral special 
conformal supersymmetry. In addition, it has conformal weight (—2) and U(l) charge (+4) 
at each point. Most importantly, these properties fix the form of X n , up to an arbitrary 
conformally invariant factor depending on the space-time variables Xi only. 

Comparing the conformal and SU(A) properties of X^ with those of the correlation 
function G®\. t defined in (I2.13p . we observe that they match except for the conformal 
weight. For the correlation function, the weight at each point is fixed by the conformal 
properties of the stress-tensor multiplet (I2.6P to be (+2) whereas for X 4+ ^ it equals (—2). 
Therefore, the nilpotent function I±+e should be multiplied by a conformally covariant 
function of Xi with weight (+4) at each point. Indeed, from ( I2.15P and (I2.16P it follows 
that for I = 1 and 1 = 1 the correlation function G^) e . t has the expected form. 

Drawing together these facts, we can write Gf^_ £ . e for arbitrary £ in the form 

G^/l, . . . , 4 + 1) = 2 ( { _ N J~H x X 4+ , x f\x u . . . , x A+i ) , (2.19) 



S 



where f^ > is a conformally covariant function of x\, . . . ,x^£ with weight (+4) at each 
point. Here the constant factor on the right-hand side of (12.1 9p is introduced to simplify the 
expression for f^' (see Eq. (I2.2ip below). Most importantly, it follows from the full crossing 
symmetry of the correlation function (12. 8p together with the aforementioned permutation 
S n symmetry of the prefactor I^+e that this function must be completely symmetric under 
exchange of any of the 4 + £ points: 



(. . . , Xi, . . . , Xj , . . .) (. . . , Xj, . . . , Xi, . . .) . (2.20) 

This is a non-trivial permutation symmetry, specific to the J\f = 4 theoryjf] We will show 
below that the S^ + e symmetry, together with the conformal and some other properties, 
impose strong constraints on the possible form of the function Our aim is to restrict 
it to a unique expression using the correlation function/amplitude duality of [3T1[32] . In 
this paper we do this in the cases £ = 3,4. 

For £ = 1 and £ = 2 we find from (12.151) and (12.161) the explicit expressions for the 
function 

1 



/ {2) (*i, 



^5 



X6) 



2 ' 

<i<j<5 X ij 



48 ^ [T 



2 .y-t 2 /y»2 
"2 "3"4 US"6 



2 

<j<j<6 X ij 



(2.21) 



These relations can be represented in a diagrammatic form as shown in Fig. [TJ 





t=\ £=2 

Figure 1: Diagrammatic representation of the function f^(x) for t = 1 and £ = 2, Eq. (|2.21|) . 
Each line with labels i and j at the end points denotes a scalar propagator 1/a^-- 

We recall that the Borm-level correlation function G^l e . e defines the integrand of the 
four-point correlation function (12.141) at £ loops. Substituting (12.191) into (12. 14j) we find 
that the £— loop correction to the four-point correlation function takes the form 

Gi £) (l,2,3,4) = 2( ( ^~ 4 1} x £(1,2,3,4) x F® for £ > 1 , (2.22) 

3 Explicit perturbative calculations [3] in a generic conformal Af = 2 theory show that in general the 
two- loop integrand does not have this symmetry, it is only there in the Af = 4 case. 
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where _R(1,2,3,4) is given by (I2.18P and the function is defined by the Euclidean 
integral 

F^(x 1 ,x 2 ,x 3 ,x,) = ^A^A^A I ^...^/W^...,^) . (2.23) 

Taking into account the conformal properties of the function we find from f!2.23j) that 
F™' is a conformally covariant function of x±, x 2 , X3, x$ with weight (+1) at each point. 
The prefactor R(l, 2, 3, 4) has the same conformal weights. Thus, their product on the 
right-hand side of ( 12 . 22 j) has conformal weights (+2) as it should be for the four-point 
correlation function G^l, 2, 3, 4). 

The presence of the universal rational prefactor R in f 1 2 . 2 2 [) for any £ signifies that the 
loop corrections to the four-point correlation function of M = 4 stress-tensor multiplets is 
determined by a single function of the four points X\, . . . ,£4. This is much less freedom 
than what one could expect in the six SU(4) channels of the tensor product 20' x 20' = 
1 + 15 + 20' + 84 + 105 + 175. This fact is known under the name of "partial non- 
renormalisation theorem" [5]. 

The relation ( I2.23P illustrates the non-trivial meaning of the permutation symmetry 
f l2.20p of the function In application to (I2.23p . it exchanges the external points 

Xi,x 2 ,x 3 ,x 4 and the integration points x 5 , . . . ,x i+ i and, therefore, puts the two sets of 
points on an equal footing. We would like to emphasise that this symmetry is obscured for 
the function F^ by the integration on the right-hand side of (12. 23 p . since it breaks the 
symmetry of the function f^> down to the S4 symmetry of the function F^K 



3 The four-point correlation function at three loops 

In this section, we combine the permutation symmetry of the integrand together with the 
conjectured correlation function/amplitude duality to work out the three-loop expression 
for the four-point correlation function (12. 9p . 



3.1 General Ansatz for the three-loop integrand 

Let us find out to what extent we can determine the rational function in (I2.23p . for 
example at three loops, for i = 3. To do this we examine again the one- and two- loop 
functions given in (12.21 j) and note that they have a very specific denominator of the form 
rii<j x lj- A s we will see in a moment, this form of the denominator can be deduced from 
the operator product expansion (OPE) analysis of the correlation function (I2.19p and it 
should be universal for any £. This suggests to write the general Ansatz for the function 
as follows: 



PW( Xl ,...,X 4+£ ) 



III 



<i<j<4+£ x ij 



where is a homogeneous polynomial in x?- which is invariant under S± + t permutations 
of Xi. We recall that the function is conformally covariant with weight (+4) at each 
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point, while the denominator on the right-hand side of ( 13. ip has weight — {£ + 3). As a 
consequence, the polynomial P^ should have uniform weight —(£—1) at each point, both 
external and internal. This implies, in particular, that P™' is a polynomial in x\j of degree 
(£ —!){£ + 4)/2. Indeed, the one- and two-loop expressions (I2.2ip have this property. At 
three loops, for i = 3, the degree of P^ should equal 7. 

A characteristic feature of the Ansatz (13. II) is that the distances between each pair of 
points appear in the denominator on the right-hand side of fl3.ll) to the first power only. 
In other words, for a generic polynomial P^> the function (13. ip scales as f^\x) ~ ^-l x % 
for Xi — >■ Xj. This property can be understood as follows. We recall that according to 
(I2.13p . (I2.17P and (I2.19P the function f^\x) determines the Born-level expression for the 
correlation function involving the protected operator and the M = 4 SYM Lagrangian 
L (itself protected, since it is a member of the same supermultiplet as 0): 

(0(1) . . . 0(4)£(5) . . . £(4 + ~ R(l, 2, 3, 4) x\ 2 x\x\Ax\A /W( x ) . (3.2) 

As a result, the asymptotic behaviour of the function /w(a?) for Xi — > Xj should follow 
from the (Born-level) OPE expansion of the operators located at the points x» and Xj. 
Depending on the choice of points, we can distinguish three different cases: 0(i)0(j), 
0{i)C{j) and £(i)£(j). 

In the first case, the OPE expansion of the product of two half-BPS operators is spec- 
ified in Eq. (14.51) below. For X{ — > Xj it is dominated by the contribution from the identity 
operator 0(1)0(2) ~ I j x\ 2 + 0(1 / x\ 2 ) ■ However, the identity operator only contributes to 
the disconnected part of the correlation function, (0(1)0(2)) (0(3)0(4)£(5) . . .£(4 + £)). 
As a consequence, the leading contribution to the connected correlation function in (13. 2p 
comes from subleading 0(l/x\ 2 ) terms in the OPE (see Eq. (14.51) ). thus implying that 
the right-hand side of (13.21) should scale as l/x\ 2 - Indeed, it follows from (I2.18P that 
i?(l,2,3,4) ~ 1/^12 and, therefore, x\ 2 f^ l \x) should stay finite for x\ — > x 2 . This leads 
to f {e) (x) ~ 1 jx\ 2 in agreement with (13.11). Let us now examine the OPE expansion of 
0(l)£(j). For x\ — > Xj the leading contribution scales as cO{l)/x\y However, the corre- 
sponding coefficient function is proportional to the three-point function c ~ (00£) which 
vanishes in M = 4 SYM to all loops (non-renormalisation of the protected two-point func- 
tion (00)). As a result, the dominant contribution to the OPE scales as \jx\^ again 
leading to the asymptotic behaviour f^(x) ~ Finally, we consider the OPE ex- 

pansion of the product of two Lagrangians. It is known that the dominant contribution 
to C(i)C(j) as Xi — > Xj comes from contact terms proportional to S^(xi — Xj) and its 
derivatives. Since the correlation function (I3.2p involves Lagrangians at distinct points, 
the contact terms do not contribute to (13. 2p and we are left with the contribution 
This immediately implies that f™\x) ~ The above considerations explain 

the form of the denominator in our Ansatz (13. ip . 

4 Thc contact terms play an important role in the Lagrangian insertion procedure based on the dif- 
ferentiation of the correlation function with respect to the coupling constant. Namely, they ensure that 
successive differentiation just inserts the action without reproducing the correlation function itself. 
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Making use of the Ansatz (13.1 j) for £ = 3, we can reduce the problem of constructing 
the rational function f( 3 '(xi, . . . , Xj) to that of finding S7 symmetric conformally covariant 
polynomial P^ 3 \xi . . . , X7) of weight (—2). As was already mentioned, takes the form 
of a homogeneous polynomial of x\j of degree 7 and it is not hard to convince oneself that 
there are in fact only four possibilities HE] for P^: 

(a) heptagon: ^i2 x 23^34 x 45^56 x 67 x 7i + $7 permutations , 

(b) 2-gon x pentagon: (^X^L^^se^CT^ra) + SV permutations , 

(c) triangle x square: { x \2 x ^ x \^{ x A5 x \§ x §7 x 7a) + S 7 permutations , 

(d) 2-gon x 2-gon x triangle: (^K^K^se^ 2 -^) + ^7 permutations. (3.3) 
It is convenient to represent the different choices in the form of diagrams as shown in Fig. [2J 
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(a) (b) ( C ) (d) 

Figure 2: Diagrammatic representation of the different S7 symmetric polynomials P^(xi) defined 
in (|3.3p . The indices (a±, . . . , 7) correspond to all permutations of the external points (1, . . . , 7). 
The dashed lines with indices cr, and aj at the end points denote the factors x\ a .. 



In the expression for the function _p 3 \ Eq. (13. II) . the polynomial P^ is divided by the 
product of distances between any pair of points. Replacing P^ in (13. ip by the various 
expressions defined in (13.31) . we find that some factors x\^ in the numerator cancel out 
against similar factors in the denominator in (13. ip . In diagrammatic terms, a solid line 
and a dashed line between two points annihilate each other. The resulting four different 
contributions to are represented graphically in Fig. [3j An important feature of these 

5 In fact, the polynomials listed in (|3.3[) are not independent in four dimensions due to a single Gram 
determinant constraint relating all four of them [33] . 

6 In principle, with six points or more one can form pseudo-scalar conformal invariants which could 
appear in the integrand. However, these involve e tensors, which are antisymmetric under interchanges 
and are thus inconsistent with the aforementioned S^+i permutation symmetry. Moreover, the four-point 
correlation function (the integral (|2.23[) itself) cannot contain parity-odd terms. This shows that, even if 
such terms exist at the level of the integrand, they must be total derivatives which do not contribute to 
the integral. For this reason we do not consider pscudo-scalars in p. 31) . 
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diagrams is that for each vertex the number of solid lines minus the number of dashed 
lines attached to it equals 4. This property ensures that the conformal weight of the 
corresponding contribution to equals 4 at each point. 



a 5 a 5 03 07 




02 <J\ 06 04 Gl °"5 



(a) (b) (c) (d) 

Figure 3: Diagrammatic representation of the four classes of functions f^\x) corresponding to 
the polynomials shown in Fig. [2j Solid lines denote scalar propagators \jx 2 a . a . while dashed lines 
stand for numerator factors x 2 „ „ . 

The general expression for the polynomial satisfying the symmetry constraints 
specified above is given by a linear combination of the four terms in ( 13. 3 p with arbitrary 
coefficients. Similarly, the general expression for p 3 > is given by the same linear combina- 
tion of four diagrams shown in Fig. [3j We would like to emphasise that this result holds 
for a gauge group SU (N c ) with arbitrary N c . 

To fix the value of these coefficients we have to impose some additional conditions. 
One possible approach is to try to extract some anomalous dimensions, for instance that 
of the Konishi operator, from the OPE analysis of the correlation function. This requires 
the ability to compute the integrals whose integrands are encoded in Fig. [3] in the singular 
short-distance limit Xi — > Xj. We explore this approach in Section H] in the simpler case of 
the planar correlation function (i.e., for large N c ), although in principle one could generalise 
it to the non-planar case as well. 

Another approach consists in exploiting the conjectured duality between scattering 
amplitudes and correlation functions [311 [32]. Since this duality only works for planar 
amplitudes and correlation functions, we are thus restricting our study to the planar sector. 
As we show in the next subsection, the comparison of the correlation function with the 
known result for the three-loop four-point amplitude allows us to unambiguously fix the 
coefficients of all four topologies in Fig. [3J 

3.2 Correlation function/scattering amplitude duality 

In application to the four-point correlation function, the duality establishes the correspon- 
dence between the four-particle scattering amplitude v4 4 in planar M = 4 SYM and the 
correlation function G4 in the limit where the four operators become light-like separated 
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in a sequential fashion, x{ 2 = ^23 = x 34 = x \i = ; 

2 lim (G 4 (x)/G 4 0) (x)) = (A 4 (p)/4 0) (p)) 2 . (3.4) 

Here on the right-hand side the amplitude Ai(pi,P2)P3)P4) depends on the light-like mo- 
menta of the scattered particles, pf = 0. They are identified with the coordinates Xi of 
the operators on the left-hand side through the dual space relation pi = xt — Xj+i with the 
periodicity condition x i+4 = X{. This duality is understood at the level of the integrands 
on both sides of the relation, and not in terms of the divergent integrals. 

We recall that the i— loop correction to the correlation function G4 is determined by 
the scalar function F^>(xi, x 2 , £3, X4) defined in Eqs. (I2.22p and (I2.23p . Similarly, the 
I— loop correction to the four-particle amplitude in the planar M = 4 SYM theory is 
given by the scalar function (pi, p 2 , P3, Pi) = Af\p)/A^\p). At present, the functions 
are known up to five loops, for £ = 1, . . . , 5, in the form of a sum of scalar £— loop 
planar Feynman integrals of various topologies [3HH301SS] - The duality (13. 4p can then be 
formulated as a relation between the integrands in and as follows. In the light-like 
limit we have, from (12.181) and (12.111) . 

p/i 9 o A\ v ™2 2 ^121/23^34^14 r<( Vl 9 Q ^ ~ 1 ^12^23^34l/l4 /o cr\ 

X 13 X 24 „ „ „ 2 ' ^4 l 1 ,^^^ /,,_2\4 ~2 „2 „2 ™2 ' K ' ) 

X 12 J '23- X '34 J '14 J x 12 x 23" ( '34 x 14 

hence 

2 lim (Gf/Gf) = lim {2x\AF^) . (3.6) 
Then the duality relation (13.41) becomes 



lim [ 1 + 2x13X34 



J^a^A = (l + J2^M (i) ) , (3.7) 

£>1 J \ £>1 J 



or equivalently 



lim (x 2 3 x 2 4 FW) = M« , (3.8) 
lim (x 2 13 xj 4 FW) = + I (M^) 2 , (3.9) 
lim (x? 3 x^ 4 F (3) ) = M (3) + M (1) M (2) , . . . . (3.10) 



c 2 



c 2 



The first two relations, Eqs. (13. 8 p and (13. 9p . have been verified in Ref. [31] by making use 
of the explicit two-loop expressions for the correlation function and scattering amplitude. 
Our task now is to do the next step and use the new relation (I3.10p to fix the coefficients 
of the four topologies in Fig. El 
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For the three-loop correlation function, we deduce from (12.231) that the integrand of 
_p(3) j s given by 

2 2 2 2 2 2 

nn" rip*- 1 ty*- 1 nn'-' ry^ nn'-' 

[F^} mtcgrand = 12 g^)3 4 X /(3)(Xl ' • • • ' X7) ' (3 - U) 

with /( 3 ) given by a linear combination of graphs shown in Fig. |3j In the light-like limit, 
xji+i 0, the prefactor on the right-hand side of (13. lip vanishes but (some of the terms 
in) the function J® develops poles in l/xf i+1 , so that their product stays finite. The 
surviving terms will appear on the right-hand side in Eq. (I3.10p . 

Let us find out which integrals we expect to see in (13.101) . The one- and two- loop 
four-gluon amplitudes, and M^ 2 \ involve one- and two-loop ladder diagrams shown 
in Fig. |U(a) and (b), respectively. The three-loop four-gluon amplitude M^> involves only 
two integral topologies, the three-loop ladder and the so-called "tennis court" [38]. They 
are depicted in Fig. H] (c) and (d), respectively, both as conventional momentum p— space 
diagrams and as dual x— space diagrams. 




(a) (b) (c) (d) 

Figure 4: Dual conformal x— integrals (solid lines) and momentum p— integrals (doted lines) in the 
three-loop planar four-gluon amplitude: (a) one-loop ladder, (b) two-loop ladder, (c) three-loop 
ladder and (d) tennis court. 

The one- and two- loop graphs are redrawn again in Fig. [5] with the external and internal 
points labelled and, in the case of the two-loop ladder, with an extra numerator factor 
(dashed line) added. This factor balances the conformal weights at points 2 and 4, so that 
the integral has uniform conformal weight (+1) at each external point. Similarly, the tennis 
court T and three-loop ladder L have been redrawn in Fig. |6] with the necessary dashed 
lines added. Fig. [6] contains another diagram of three-loop topology, the product g x h of 
the one- and two- loop ladders that we expect to find in the non-linear term M^M^ on 
the right-hand side of (13. IQj) . 

We may interpret the graphs in Fig. El with conformal weight (+1) at each external 
point, as those terms in the function F^ 3 \ Eq. ( 13. lip , which survive in the light-like limit. 
Now, we are interested in the manifestly SV symmetric function /( 3 ). In order to up- 
grade the F— terms shown in Fig. E]to /—terms, we need to divide them by the prefactor 
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0(1,2,3,4) 



3; 2, 4) 



Figure 5: Diagrammatic representation of the one- and two- loop integrals (|3.13p . The solid lines 
denote scalar propagators and the dashed line denotes a factor of in the numerator. Each 
external point has conformal weight (+1). 




T(l,3;2,4) L(l,3;2,4) g x h(l, 3; 2, 4) 

Figure 6: Diagrammatic representation for the three-loop integrals (|3.14p surviving in the light- 
like limit. The dashed lines (numerator factors) ensure uniform conformal weight (+1) at each 
external point, as well as at the integration point 7 for T. 



x i2 x i3 x i4 x 23 x 24 a '34 from ( 13 .11)) . In graphical terms, this means to superpose a square with 
diagonals (all lines solid) and with vertices at points 1,2,3,4, onto each graph in Fig. El In 
the process some numerator factors (dashed lines) cancel against propagators (solid lines). 
The result of this manipulation are graphs of the type shown in Fig. |3j It is not difficult 
to see that all three F— graphs from Fig. [61 when upgraded to /—graphs, fall into the 
topology Mh). 

We conclude that the three other topologies shown in Fig.[3Ja),(c),(d) should not appear 
in the light-like limit, otherwise they would give unwanted contributions to the amplitude. 
It is not hard to check that none of the topologies in Fig. |3] entirely vanishes in the light-like 
limit. Hence, the amplitude/correlation function duality leads us to setting the coefficients 
of the topologies in Fig. [3Ja),(c),(d) to zero0 The only coefficient which still needs to be 
fixed is that of the topology in Fig. [3Jb). We do this in Section 13.41 by examining the 
details of the duality relation (13. 10)) . 

We wish to point out that there exists an alternative way to interpret the duality relation 

We repeat that this only concerns the correlation function in the planar limit, to which the duality 
applies. In the general case (for an arbitrary gauge group SU(N C )) all four topologies in Fig.|3]are expected 
to contribute to the correlation function. 
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(13. 7p . We may say that the entire planar amplitude or correlation function is generated, 
up to three loops, by the lower- loop non-linear terms (M^) 2 in ( 13. 9 p and M^M^ in 
( 13.1 Op . Indeed, at one loop the only possible choice for the numerator of f^ 1 ' in ( I2.2ip is 
a constant whose value is fixed by the relation (13. 8p . Then, at two loops the non-linear 
term (M^) 2 in appears, together with (with exactly these coefficients), in the 
unique symmetric choice for in (12.211) . This means that the non-linear term actually 
tells us what the linear two- loop term must be, and hence determines the two-loop 
amplitude. Similarly, the non-linear term M^M^ in (13.101) determines the three-loop 
term M^ 3 \ by the simple fact that they belong to the unique symmetric topology 3(b) 
contributing to f^ 3 \ In Sect. [6] we comment further on the extent to which higher-loop 
amplitudes/correlation functions can be determined from lower-loop amplitudes in this 
way more generally. 



3.3 Three-loop integrand 

As was explained in the previous subsection, the three-loop integrand for the planar 
four-point correlation function can only receive contributions from the topology shown 
in Fig. E](b). The corresponding expression for the function /( 3 ) reads (up to an overall 
normalisation factor) 

1 4 2 2 2 2 2 



III 



2 

<i<j<7 X ij 



where the factor 1/20 simply ensures a sum over distinct terms: there are 7!/20 = 252 
distinct terms in the sum over Sj permutations. 

Defining the three-loop integrand (13.111) we have to split the seven points Xi into four 
external points Xi,...,X4 and three internal (integration) points x^,xq,X7. Then, the 
SV symmetric rational function J® can be decomposed into a sum of terms separately 
invariant under S4 x S3. Here S4 and S3 act on the external and internal points, respectively. 
This gives the five distinct terms displayed in Figure [7J Notice that, drawn like this, all 



(a) (b) (c) (d) (e) 

Figure 7: Diagrams contributing to the function /( 3 ) in the three-loop integrand. White nodes 
represent the external points x±, . . . , X4 and black nodes represent the internal points £5, xq, £7. 



17 



the graphs in Fig. [7] have the same topology, the only difference being the localisation of 
the external (white) and internal (black) points. 

Recall that the integrand (13. lip includes an additional factor given by the product of 
distances between all external points. This factor breaks the S7 symmetry of the function 
/( 3 ) down to S 4 x S 3 . At the same time, it simplifies the form of the graphs contributing 
to by removing the extra propagators between the external points. Some of the 
corresponding integrals, those which survive the light-like limit (for x\ 2 = x\\ = x\± = 
%1i — 0), already appeared in Fig. O Away from this limit, for arbitrary xf i+1 7^ 0, we find 
only two new conformal three-loop integrals E and H (for "easy" and "hard" , referring to 
their Mellin-Barnes evaluation in Appendix [Cj , which are depicted in Fig. [HJ Thus, the 
complete list of integrals which we expect to find in the correlation function (away from 
the light-like limit) is shown in Figs. E] and [SJ All others graphs for are obtained by 
permuting the external points. 




£(1,3; 2, 4) #(1,2; 3, 4) 

Figure 8: Diagrammatic representation for the new three- loop integrals ()3.14[) contributing to the 
four-point correlation function away from the light-like limit. The dashed lines between adjacent 
external points indicate that these integrals vanish in the limit. 

Here are the expressions for the one- and two-loop ladder (box and double-box) integrals 
g and h from Fig. 



0(1,2,3,4) 



d A x§ 



A /7T-2 / /"V«2 /y»2 ,-y»2 , y» 2 ' 

<±/l ./ X 15 X 2 5X 35 X 45 



h(l 2- 3 4) = X ™ [ d 4 x 6 

^ ' ' °' > (Att2\2 J (,,.2 „2 ™2 \™2 1^2 J J \ ' yo.±oj 
\f±li ) J yx 15^35^45 M56 v x 26 x 36 J ' 46/ 
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and for the three-loop integrals from Figs. M and |H1 



T(l,2;3,4) 
£7(1,2; 3, 4) 
£(1,2; 3,4) 
(<7X /*)(!, 2; 3, 4) 



x 34 f d^x^XQcl^Xj x\ 7 



A '7T-2 ^3 / ( - y» 2 ,-y» 2 \ /^y,2 ,->-» 2 \ / 2 /y» 2 ,->-» 2 \ ™2 ™2 ,-y»2 

> 4,1 J i 1^15^35 A x 16 x 46A x 37 x 27 x 47J- z, 56 J/ 57-' f '67 



_/l -77-2^3 / ry>2 ~, 2 ™ 2 \ ^,2 / ™2 ^2 ™2 \ ,-vj2 / ™,2 ™2 ™2 \ ' 

^71 ) j yx 15 x 25 x 35 )x 56 {x 26 x 36 x 46 )x 67 {x 17 x 27 x 47 ) 
xi /" d 4 x 5 d 4 x G d 4 xj 



/|'77-2^3 / ( ry 2 /y» 2 ™ 2 \ ,-y»2 / ™2 ^v>2 \™2 / ^»2 ™ 2 ™ 2 "\ 
^71 j J ^15^35^45^56 l x 36 x 46>' X 67V X 27 X 37 X 47 ) 

12 x 34 



X 12 X 34 /" c? 4 x 5 (i 4 x 6 (i 4 x 7 



/I -77- 2 ^3 / ( i~y* 2 , >-< 2 ,y»2 ^v«2 \ / ™2 rp2 ,-y»2 \ / ,-y. 2 ,-y»2 . _ \ ,.2 
^ n ) J l x 15 x 25 x 35 x 45A x 16 x 36 x 46A x 27 x 37"''47M67 
„2 „.2 ^2 /■ J4^ j4„ J4^. „.2 



£Tfl 2-3 4) x 4:i x 23 x 34 f d x^d x^d x-r x 57 . 

^ /| ^7-2^ 3 / ( |-V»2 j-»2 y>2 /y»2 \ ,7^2 / /y»2 ,-y»2 \ ,-y. 2 / ,y>2 /y>2 ,-y» 2 ™ 2 \ V 

^ ) j \,x 15 x 25 x 35 x i5 )x 56 yx 36 x 46 )x 67 {x 17 x 2 7X 37 x /k7 ) 

Summarising our analysis, we obtain the full four-point correlation function up to three 
loops as follows: 

G 4 (l, 2, 3, 4) = Gf + 2 ^ ~ 4 1} 2, 3, 4) [aF« + a 2 F^ + a 3 £^ + 0(a 4 )] , (3.15) 

where F« and can be read off from fF2T2T|) . and fj57T5j) . 

= (7(1,2,3,4), 

F (2) = /i(l, 2; 3, 4) + h(3, 4; 1, 2) + /i(2, 3; 1, 4) + h(l, 4; 2, 3) 

+ /i(l, 3; 2, 4) + fc(2, 4; 1, 3) + | {x\ 2 x 2 u + s? 3 a& + x 2 u x 2 23 ) [g{l, 2, 3, 4)] 2 , (3.16) 

and F( 3 ) is given by the sum over all topologies shown in Figs. |6] and [HJ symmetrised over 
all 5*4 x S3 permutations of white and black nodes. The explicit expression for reads 

F (3) = ^(l, 3; 2, 4) + 11 perms] + [£7(1, 3; 2, 4) + 11 perms 

+ [£(1,3; 2, 4) + 5 perms] + [(g x h)(l, 3; 2, 4) + 5 perms] 

+±[#(1,3; 2, 4) + 11 perms], (3.17) 

where "+11 perms" etc. denotes a sum over the remaining distinct 5*4 permutations of 
external points 1,2,3,4. The reason that the hard integral H comes with a factor of a 
half is simply that its integrand has the additional symmetry under £5 -H- £7 and we 
are summing over all distinct integrands with weight 1 (and then dividing by 3!). One 
can easily check that the total number of integrands on the right-hand side of (I3.17P is 
3! x (12 + 12 + 6 + 6 + 12/2) = 252, as it should from (I3TT2"]) . We recall that the expression 
for F^ is defined up to an overall coefficient. As we show in Section 1374"} its value is fixed 
from the duality correlation function/amplitude to be 1. 
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There exist identities that certain conformal off-shell integrals (but not their inte- 
grands!) satisfy for xf i+1 ^ (see Ref. [IS])- In Section |4.2. II we will use these identities to 
simplify the expression (13.171) further. It should be noted, however, that all of the above 
integrals diverge in the light-cone limit xf i+l — > and need regularisation. In this limit 
the conformal identities do not hold anymore, so for the purpose of checking the duality 
correlation function/amplitude we ought to treat all integrals in (13.171) as independent. 

3.4 The square light-like limit and the four-gluon amplitude 

Let us now verify that the proposed expression for the correlation function (13. 17ft reproduces 
the known result for the three-loop amplitude in accordance with the amplitude/correlation 
function duality (I3.10p . In the process we will also fix the overall normalisation (the 
coefficient of topology 3(b)). 

To this end, we go in Eqs. (13 . 1 5[) - (13. 1 7[) to the light-like limit x\ 2 , x% 3 , x| 4 , x\ x — > and 
observe that the expression on the right-hand side of (13.171) reduces dramatically. As was 
mentioned in the previous subsection, all permutations of the integrals E and H vanish 
in this limit and the only non-vanishing contributions come from 4 permutations of T, 2 
permutations L and 2 permutations of g x h. Explicitly we obtain 

lim =0(1,2,3,4), 

lim F& = h(l, 3; 2, 4) + h(2, 4; 1, 3) + \x\ 3 x\± [g(l, 2, 3, 4)] 2 , 

lim F< 3) = T(l, 3; 2, 4) + T(l, 3; 4, 2) + T(2, 4; 1, 3) + T(2, 4; 3, 1) 

+ L(l, 3; 2, 4) + L{2, 4; 1, 3) + ( 3 x h){l, 3; 2, 4) + ( 9 x h)(2, 4; 1, 3) . (3.18) 

We recall that the integrals in (13. 18[) are divergent and require regularisation. However, for 
the purpose of testing the duality relation (I3.10p we will only need their integrands which 
are well-defined rational functions of the coordinates of the four external points and the 
three integration points. 

Let us now compare (I3.18P with the analogous expressions for the four-gluon amplitude. 
The latter are naturally defined as functions of the on-shell particle momenta pi. Going 
to the dual space representation, pi = Xi — x^+i, the one-, two- and three-loop four-gluon 
amplitudes are given by [16] 

M[ 1] = 4,4^(1,2,3,4), 

Mj 2) = xj 3 x 2 M [h(l, 3; 2, 4) + h(2, 4; 1, 3)] , 

Mf = x\A [T(h 3; 2, 4) + T(l, 3; 4, 2) + T(2, 4; 1, 3) + T(2, 4; 3, 1) 

+ L(1,3;2,4) + L(2,4;1,3)], (3.19) 
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with the integrals defined in (13.131) . (I3.14p . Inserting (I3.19P into (13.101) and comparing with 
f)3.18p we find precise agreement with the amplitude/correlation function prediction, as 
stated. 

This completes our construction of the three-loop four-point correlation function based 
on the usage of the hidden SV symmetry of the integrand combined with the correlation 
function/amplitude duality. 

4 OPE test of the three-loop correlation function 

In this section we check that the proposed three-loop four-point correlation function ( 13 . 1 7[) 
agrees with the operator product expansion (OPE). For protected scalar operators of the 
type (12.1 1) it takes the following general form (here we do not display the SU(4) index 
structure): 



0(x 1 )0(x 2 ) = \Cooo A —r-^ UX ^-Ax 12 U • • • (x 12 ), s [0^ s {x 2 ) + ...]. (4.1) 




The sum on the right-hand side runs over conformal primary operators O^ 1 "'^ carrying 
Lorentz spin S and scaling dimension A and over their conformal descendants shown by 
dots. Here we took into account that the scaling dimension of the operators (12.41) equals 
2 and it is protected from quantum corrections. The contribution of each operator to the 
right-hand side of (14. ip is accompanied by the coefficient function Cooo& ■> which determines 
the three-point correlation function (0(xx)0(x 2 )O l £"' tJ ' s (x 3 )) . 

In this section, we shall apply the OPE (14.11) to reproduce the known result for three- 
loop anomalous dimension of the Konishi operator (see Eq. (14 .2p below) from the four- 
point correlation function (13.171) and to predict its three-point correlation function with 
the protected scalar operators (12.41) . 

4.1 OPE expectation 

Let us examine the asymptotic behaviour of the four-point function G 4 (l,2,3,4) in the 
double short- distance limit X\ — > x 2 , x% — > £4. We would like to emphasise that in this 
section we consider the correlation function in the Euclidean regime, so that the limit 
x\ 2 —> is indeed a short- distance (coincidence) limit, Xi — > x 2 , and not a light-cone limit 
like in Section [3.41 The advantage of the double short- distance limit, x\ 2 ,x\ A — > 0, is 
that we can apply the OPE simultaneously to two pairs of operators in G^l, 2,3, 4) and 
relate the leading asymptotics of the four-point correlation function to the contribution of 
a particular conformal primary operator, the Konishi operator. 

It follows from (14. ip that for x\ 2 ~^ the leading contribution to the OPE comes 
from the operators with the minimal twist = (A — S) and the minimal Lorentz spin S. 
In particular, the identity operator X (with A = S = 0) produces the most singular 
contribution ~ l/xf 2 . The first subleading 0{l/x\ 2 ) contribution originates from the 
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gauge-invariant operators carrying zero Lorentz spin 5 = and tree-level (naive) scaling 
dimension A^ ) = 2. In J\f = 4 SYM such operators are built from the six scalar s $ 7 . They 
have the general form C/jtr($ 7 $ J ) with some coefficients Cu and carry 5£7(4) irreps from 
the tensor product 6x6 = 1 + 15 + 20'. In this way, we identify two such operators: 
the half-BPS operator O^q,, Eq. (12. ip . belonging to the 20' (or the symmetric traceless 
rank-2 representation of 5*0(6)) and the singlet Konishi operator^] 

JC = tr($ 7 $ 7 ) . (4.2) 

The scaling dimension of the half-BPS operator O 2 o' is protected from quantum corrections 
and takes its canonical value A@ = 2, while the Konishi operator acquires anomalous 
dimension: 

oo 

Ak = 2 + 7 ,c(a) = 2 + ^aV (4-3) 

i=i 

The goal of our test is to extract this anomalous dimension from the four-point correlation 
function up to three loops and to compare it to the known values from the literature 

nnmsz]: 
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7i = 3 , 72 = -3 , 73 = — . (4.4) 

Thus, restricting to the contributions of operators with naive scaling dimension up to 
two, the OPE of two operators 0(x,y), Eq. (12. 4p . contains the operators X, K and O 2 o r - 

0(x u yi) 0(x 2 , y 2 ) = ci ^ Yl J 2 ^ X + c K (a) ^ JC(x 2 ) 

+ co^^-YuYij 0&(x 2 ) + . . • , (4.5) 



x 



12 



where Y\ and Y 2 denote the 5*0(6) harmonics depending on y 1 and y 2 , respectively (see 
Eq. (lA.ip in Appendix EJ. The polynomial harmonic structure on the right-hand side of 
( 14. 5 p is determined by the usual requirement of harmonic analyticity and by the SU(4) 
weights (+4) at points 1 and 2 (each Y carries weight (+2)). The dots in (14. 5 j) denote 
the contribution of operators of higher spin and twist. It is suppressed by powers of x± 2 
in the (Euclidean) limit x\ — > x 2 . The constants cj and cq determine the two-point and 
three-point functions of the half-BPS operators 02O' ; which are known to be protected. 
Therefore, the constants c% and cq do not depend on the coupling constant and are given 
by their tree-level expressions 

A^ 2 — 1 1 



8 Thc irrcp 15 in the tensor product, i.e. the antisymmetric rank-2 tensor (the adjoint of SU(4)), cannot 
be realised as a scalar bilinear operator since tr($ 7 <I> J ) is symmetric in / and J. 
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In the double short-distance limit X\ — > x 2 , x% — > £4 we can apply the OPE (14. 5 p to 
the products of operators 0(xi, yi)0(x2, 2/2) an d 0(x^, yz)0{x* i Va) an d express their four- 
point correlation function in terms of the two-point functions | 

<0&(* 2 ) 0%( Xl )) = (s' K S' L + S'V - \6"6 KL ) . (4.7) 



^24 

Unlike the half-BPS operator O^i, the Konishi operator K. is not protected, so the nor- 
malisation constants die and the coefficient c^(a) on the right-hand side of (14. 5 p depend on 
the choice of the renormalisation scheme. We shall adopt the convention that the Konishi 
operator is normalised such that d/c keeps its free-field value: 

^ = 3^, c K (a) = T ±- 2+ 0(a). (4.8) 

Below we show (see Eqs. (I4.13P and (14.150 ) that the OPE allows us to find c^(a). 

Taking into account these relations and using ( 14. 5 j) we finally obtain the asymptotic 
behaviour of the four-point correlation function for x\ — > X2, £3 — > £4 as 

4 4 

n _ VyiVza 2 

( - T 4 — —7 T~ Cz 

x 12 x u 

J '12 J/ 34 J '24 z 

+ j¥^r (4W«fc «™ w/2 - I <&*) + ... , (4.9) 

J '12" t '34 x 24 V u / 

where it = ^12^34/(^13^24) is a conformal cross-ratio and the dots denote subleading terms. 
We notice using ( 14. 6 p and ( 14. 8 p that the expression in the parentheses in the last line of 
(14. 9 p scales as 0(a). 

Let us compare the OPE result (14.91) with the asymptotic behaviour of the obtained 
expression for the four-point correlation function 

G A = Gf + a G ( l ] + a 2 Gf } + a 3 Gf ] + 0(a 4 ) . (4.10) 

Replacing the tree-level expression by (12. lip and the loop corrections by (I2.22p we find in 



9 To verily the second relation in (|4.7[) , it suffices to project (|4.7p with harmonics and to take the vev 
of both sides of (|4.5p . 
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the limit x% — > x 2 , x 3 — > x 4 
G 4 = 
+ 

+ 



(^ - 1) ^ 

4(4tt 2 ) 4 x 4 2 xf 4 

(47T 2 ) 4 x l2 x 3i X 2i 

2(iV 2 -l) ^ 



(4vr 2 ) 4 * 2 2 *| 4 



[^a^w] + ... . (4.11) 



Arriving at this relation, we took into account that the prefactor R, Eq. (12.1 8p . scales at 
short distances as R(l, 2, 3, 4) = y\ 2 y\J (^12^34) + • • • • 

Comparing (14. lip with ( 14. 9 p we find with the help of ( 14. 6 p and ( 14. 8 p that the first two 
lines in the two relations coincide. Matching the expressions in the last line we obtain the 
following relation 

2(iV 2 -l)^ £ ^ 1 fo, „. 7r (a)/2 1 



(4?r 



2^4 



11 £ a^\x) = -i- U{a)d K - ~ clc x ) + . . . , (4.12) 

01 X24 V D / 



where the dots denote terms vanishing for x\ — > X2, £3 — >■ x 4 . 

Expanding the right-hand side of (I4.12p in the powers of the coupling constant and 
matching the coefficients of a e , we can obtain the OPE prediction for the short-distance 
behaviour of the functions F lyt \x). It is convenient to write 

( 1 \ N 2 — 1 

4(aMc=( 5 + t^t- (4.13) 



01 



Then, we find from (14.121) the following relations for the functions 



x 1 



x 24 F (1) (x) = —7i In u + -an 



1 

48 (i v ~" "' 



xt 4 F^(x) = — 7l 2 (Inn) 2 + — 72 + 7 7i«i hiu+ -a 2 + 



x* M F®(x) = ^7? (Inn) 3 + (i_ 7l72 + l 7l 2 onJ (Inn) 2 

11 1 \ 1 . . 

^73 + ^72"i + ^7i"2 I In u + -a 3 + . . . , (4.14) 

where the coefficients 7 ^ define the perturbative corrections to the anomalous dimension 
of the Konishi operator, Eqs. ( 14. 3 p and ( I4.4p . Notice that in the expansions on the right- 
hand sides of (I4.14p we have neglected the powers of the conformal cross-ratio u because 
it vanishes for x\ — > x 2 and £3 — > x 4 . 
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The relations f!4.14p provide powerful constraints on the form of the loop corrections. 
In the next subsection we work out the asymptotic expansion of the obtained expression 
for the three- loop four-point correlation function for x\ — > X2, £3 — > £4 and demonstrate 
that it is in agreement with ( I4.14p . Moreover, matching the coefficients of the powers of 
In u, we reproduce the known result for three-loop anomalous dimension of the Konishi 
operator (14.41) and we obtain the following values for the coefficients cm defined in (I4.13P 1^1 

3 7 /25 \ 

ai = -l, a 2 = -C(3) + -, a 3 = -f T C(5) + 2C(3) + 12j . (4.15) 

Notice that our reconstruction of the three-loop correlation function was limited to its 
leading large N c part and, therefore, the constants at in (I4.15P can in principle receive 
corrections in 1/N% starting from three loops. It is worth mentioning however that the 
Konishi anomalous dimension 7^ is known to be N c — exact up to three loops jUHTJHT]. 

As an interesting byproduct of our study of the four-point correlation function we can 
predict the three-loop result for the correlation function of two half-BPS operators O and 
the Konishi operator /C. Indeed, it follows from the OPE, Eqs. ( 14. 5 p and ( 14 .7p . that this 
correlation function looks as 

(Ofa,»0Ofa, *)*(»>)> = ■ (*■!«) 

Here the constant C(a) = cjc(a,)djc takes the form 

C(a)=fl + 3|>aM 7 ^ r , (4.17) 

with the parameters ai given up to three loops by (I4.15p . We recall that the Konishi 
operator is defined so that it keeps its free- field normalisation, Eqs. (14. 7p and (14. 8p . 



4.2 Double short-distance limit of the correlation function 

In this subsection we work out the asymptotic expansion of the three-loop four-point 
correlation function G4 in the limit xi — )■ X2, X3 — > X4. According to (I2.22p . the loop 
corrections to G4 are controlled by the functions which are given by the sum of scalar 
t— loop Feynman integrals of different topology. We start by simplifying the expression 
for the three-loop integrals defined in (I3.14p . (I3.13P and then proceed to their numerical 
evaluation using the Mellin-Barnes method. 



4.2.1 Simplifying the correlation function 

There are various identities satisfied by the off-shell (that is, for xj i+1 7^ 0) conformal 
integrals [46] allowing the expression ( I3.17P to be immediately simplified (at the level of the 

10 We acknowledge the help of Volodya Smirnov in evaluating a.3. 
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off-shell integral only, not the integrand!). They follow from the fact that the conformally 
covariant integrals depending on four external points are functions of two conformally 
invariant variables, the cross-ratios 

2 2 2 2 

u = ^ff, v = ^ff. (4.18) 

x 13 x 24 x 13 x 24 

For example, we can rewrite the one-, two- and three-loop ladder integrals g, h and L from 
(EH, flXT3j) as follows: 

X 13 X 2 4 

Ml,2;3,4) = - f i^$( 2 )(n, W ), 
x 13 x 24 

L(l,2;3,4) = ^^& 3) (u,v). (4.19) 

x 13 x 24 

All of these are finite four- dimensional integrals having conformal weight one at each outer 
point. This weight is carried by the prefactors 1/(^13^24) 011 the right-hand sides in (14.191) . 
which makes the functions ^ n \u^v) conformally invariant 1^1 Similarly, the three-loop 
integrals T, E and if from (I3.14p can be written as 

T(l,2;3,4) = ^i^T( 3 )(n, W ), 

X 13 X 2 4 

x 13 x 24 

if (1,2; 3, 4) = -Jl_H^ 3 \u,v) . (4.20) 

x 13 x 24 

Notice that the prefactor 1 / (a:^ 3 a; 24 ) and the cross-ratios u, v are invariant under simulta- 
neous exchange of the external points 1 -H- 3 and 2 -H- 4. Therefore, the integrals above 
must also be invariant under this transformation. This property is not immediately obvious 
since it is not a symmetry of the integrand itself, and it has been used in to show 
that the two-loop ladder h satisfies the "flip identity" 

fc(l,2;3,4) = fe(3,4;l,2) , (4.21) 

in addition to the manifest invariance, at the level of the integrand, under the exchange 
1 <H- 2 and 3 <H- 4. The three-loop ladder L and the "hard" integral if, in addition to the 
manifest symmetries of the integrand, satisfy the same type of flip identity: 

L(l, 2; 3, 4) = L(3, 4; 1, 2) , if (1, 2; 3, 4) = if (3, 4; 1, 2) . (4.22) 



n The ladder functions <£ >< - n ) (u, u) are known to all orders and are expressed in terms of polylogs of 
maximal degree 2n [48] . 

12 Except for the one-loop ladder function g(l, 2, 3, 4) which is manifestly totally symmetric under all 
interchanges of xi, . . . , X4 at the level of the integrand. 
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The "easy" three-loop integral £ has two non-obvious symmetries: the manifest integrand 
3 <H> 4 symmetry together with the flip identity implies 

£(1,2; 3, 4) = £(2,1; 3, 4). (4.23) 

The point labels thus come in pairs as with the other integrals so that 

£(1,2; 3, 4) = £(3, 4; 1,2). (4.24) 

Finally, the property of the two-loop ladder (I4.2ip has been applied in j4B] to flip a two-loop 
subintegral in either the three- loop ladder or the tennis court, thus proving that they are 
identical off shell (that is, for xf i+1 ^ 0) 

T(1,2;3,4) = L(1,2;3,4). (4.25) 

The identities listed above imply relations between various integrals and allow us to 
simplify the three-loop expression for the correlation function (I3.17P leading to 

£( 3 ) = 2(7(1,2,3,4) [x? 2 :4 Ml, 2; 3, 4) + /i(l, 3; 2, 4) + /i(l, 4; 2, 3)] 

+ 6 [L(l, 2; 3, 4) + L(l, 3; 2, 4) + L(l, 4; 2, 3)] 

+ 4 [£(1, 2; 3, 4) + £(1, 3; 2, 4) + £(1, 4; 2, 3)] 

+ (1 + l/v)H(l, 2; 3, 4) + (1 + u/v)H{\, 3; 2, 4) + (1 + u)H{\, 4; 2, 3) . (4.26) 

4.2.2 Mellin-Barnes representation of conformal four-point integrals 

As mentioned earlier, the one-, two- and three-loop ladder integrals g, h and L (as well 
as the tennis court integral T, see (I4.25P ) are known explicitly as functions of the two 
conformal cross-ratios [IB]- However, the new integrals £ and H we encounter in our 
analysis of the three-loop correlation function have not been studied in the literature. It is 
beyond the scope of this paper to try to evaluate these integrals analytically. Instead, we are 
going to compute them numerically in the singular double short- distance limit described 
above, in order to exhibit the relevant logarithmic behaviour expected from (14.141) . In 
terms of the conformal ratios (I4.18p . the double short-distance limit X\ — > X2, X3 — > x^ 
translates into 

u — Y , v ->■ 1 . (4.27) 

We apply the standard Mellin-Barnes method [H] adapted to the conformal integrals at 
hand. This assumes using dual x— space variables instead of the more familiar momentum 
space variables in the literature on Feynman integrals for amplitudes. 

We illustrate the procedure with the simplest example of the one-loop ladder function 
g(l, 2, 3, 4), Eq. f)3.13p . Making use of the conformal invariance in x space, we may eliminate 
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one of the propagators in the integral for #(1,2,3,4) by sending one point to infinity. For 
example, for x 4 — > oo we get from (13.131) 



hm ^g(l,2,3,4) = ~ / = ^W#,#) • 

x^oo 4vr 2 J x( 5 x 2 25 xi 5 xj 3 \xj 3 x( 3 



To determine the function of conformal cross-ratio v) from this relation, it is suf- 

ficient to work out the MB representation for the integral in (14.281) and then to replace 
x\ 2 — > ux\ 3 and x\\ —> vx\ 3 . We find (the details can be found in Appendix ICl) 

1 r io ° dz d? 

$ {1 \u,v) = -- / -^f [r(-z l n-z 2 ni + z l + z 2 )} 2 u^\ (4.29) 

where the integration z\— and z 2 — contours run along the imaginary axis and satisfy the 
conditions Rezi, Rez2 < and Re(l + z\ + z 2 ) > 0. 

The integrations in (I4.29P can be performed analytically leading to the well-known 
result [IB] for $W. For our purposes, however, we need the leading asymptotic behaviour 
of the MB integral (I4.29P in the limit ( I4.27p . It can be easily obtained from (I4.29P by 
closing the z\— integration contour to the right half-plane and picking up the coefficient of 
the double pole at z\ — 0. In this way, we arrive at 

$( 1 )(u,l) = -lnu-- + 0(u). (4.30) 



The same technique can be applied to the remaining two- and three-loop integrals in (I4.26p . 
We have summarised some steps of our evaluation in Appendix O 

Replacing the integrals in 04.26)) by their asymptotic expansions ( IC.6P we finally obtain 
the following result for the one-, two- and three-loop correlation functions: 



.4 



1 . 1 



X - 4 F (1) = - In it - - + 0(u) , 



x 4 
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X 



4 2 

F m = A ( i nu) 2_ lnu + | c(3) + 7 + 0(u)j 

F(3) = l (lnM)3 "i (lnM)2+ G c(3) + S) lnM 

^C(5) + C(3) + 6^1 + 0{u). (4.31) 



We notice that these expressions are in perfect agreement with the OPE prediction (14.141) . 
They allow us to reproduce the well-known result ( 14.41) for the Konishi anomalous dimen- 
sion and to make a prediction (14.151) for the three-loop normalisation coefficients ai defined 
in (Q5j) . 
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5 Correlation functions at higher loops 



In the previous sections, we have used the permutation symmetry of the integrand combined 
with the correlation function/amplitude duality to construct the three-loop four-point cor- 
relation function. It is fairly straightforward to extend our analysis to higher loops. As an 
example, we consider in this section the four-loop four-point correlation function. 

From (12.191) we can write the four-loop integrand for the four-point correlation function 

• • • , 8) = 2 ( ( ^ 2 ~ 1} x X 8 x /Wfo, . . . , x 8 ) , (5.1) 

where f^(xi, . . . ,x$) is conformally covariant with weight (+4) at each point and is S$ — 
invariant. According to (13.11) . this function has the following form 



P ( >(xi, . 

P'{xi...,x 8 ) = — -2-, (5.2) 

lll<i<j'<8 x ij 



where P^ is a homogeneous polynomial in • of degree 12, having conformal weight (—3) 
at each point. 

Let us analyse the possible choices for the polynomial P^(x\, . . . , x§). As at three 
loops, it is convenient to use a diagrammatic representation of P^ as a graph with vertices 
corresponding to the points x\ and the dashed lines denoting the terms x\y The fact 
that P <yA \x\ 1 . . . , Xs) has weight (—3) at each point implies that there are exactly three 
lines attached to every vertex. In graph theory terminology, such graphs are known as 
regular multi-graphs of degree 3. o Thus the problem of analysing all possibilities for 
the polynomial P( 4 \xi, . . . ,x§) is equivalent to the mathematical problem of counting all 
non-isomorphic (connected and disconnected) multi-graphs with 8 vertices of degree 3. By 
doing this, with the help of (50] , we find the total number of independent possibilities for 
P^ is 32. 

Thus, using the general properties of (15.1)) alone, we find that the integrand of the 
four-loop correlation function is parametrised by 32 unknown coefficients (one per each 
inequivalent contributing topology). To determine these coefficients we have to impose 
additional conditions on the function f^(xi, . . . , xg). So far our considerations were valid 
in M = 4 SYM with gauge group SU(N C ) for any number of colours N c . To formulate the 
above conditions, we plan to use the correlation function/amplitude duality, which is valid 
in the planar limit only. Therefore we shall restrict ourselves to the planar limit of the 
correlation function. As at three loops, the coefficients in front of the 32 topologies can 
be determined from the correlation function/amplitude duality, together with the known 
four- loop four-gluon amplitude result [39] . 



13 In graph theory a "graph" should have only 1 line between any two points, for multi-graphs we are 
allowed multiple lines between two points. Here we have such topologies, they simply correspond to 
allowing x\; in the numerator. 
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Comparison with the four-particle planar amplitude uniquely determines the four-loop 
integrand of the four-point correlation function to be 



where the contributions from topologies A, B and C are given by 

t-> I \ ■'"2222222226 

r A\X\, ■ ■ ■ ,^8) — 2^ X 12 X 13 J '16- Z '23 X 25 X 34 X 45 X 46' Z '56' Z '78 > 

p / \_^2222222244 
r B\X\ ■ ■ ■ ,Xs) — ^X 12 X 13 X 16 X 24: X2^X 3A X 38 X 45 X 5 qX 78 , 

p / \ ^222222222222 /r a\ 

C \X\ ■ ■ ■ ,X$) — Y^ X 12 a '15 a; 18 a; 23 a; 26 a '34 a; 37 a '45 a; 48 a; 56 X 67 X 78 ■ l°- 4 J 

Here the multiplicative fraction on the right-hand side are inserted to ensure that each 
term is only counted once in the sum over permutations in (15. 3p . The fact that Pa, Pb 
and Pq come with equal weight and Pq with a minus sign follows from comparing with the 
known result for the four-loop four-gluon amplitude [39] where all integrals appear with 
equal or opposite weight. 
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Figure 9: Diagrammatic representation of the three conformally covariant polynomials defined in 
(|5,4p . These are the only Sg topologies occurring in the four-loop four-point correlation function 
in the planar limit. 



The three S% topologies of the polynomial P^ which appear on the right-hand side of 
( 15. 3 P are shown graphically in Figure [9j The corresponding contributions to the function 
Eq. ( 15. 2p . are shown in Figure [TTJI 

As at three loops, all the standard integrand graphs known from the four-point ampli- 
tude are obtained from these by choosing 4 of the 8 vertices to be external vertices, and 
deleting all edges between these vertices (corresponding to external propagators). 

We would like to emphasise that the new permutation symmetry of the integrand 
function explains in part why the various four-loop integral topologies appearing in 
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Ia Ib fc 

Figure 10: Diagrammatic representation for the function f^(xi), Eq. (|5.2|) . corresponding to 
the polynomials depicted in Fig. [9l To simplify the figures we did not show the dashed lines 
corresponding to xt 7 a 8 f° r Ja and to x^ r5(Tf ,x'^ 7IJg for fs- To obtain the four-loop integrand of 
F^\xi) from these graphs, one should separate the eight vertices into four external and four 
internal points and remove the propagators connecting the external points only. 



the amplitude come with coefficients (±1). The total of 8 such integrals from Ref. [39] 
fall into the three classes A, 5, C described above. Then, the symmetry of implies 
that all the integrals from the same class must come with equal coefficients. However, this 
symmetry does not explain the relative coefficients between the three classes (for more 
comments see Sect. [6]). 

The only additional thing we need to check is that none of the 32 topologies completely 
vanish in the light-like limit (in which case we could clearly not determine them by this 
method). To prove this it is enough to show that for any allowed numerator topology 
we can choose vertices in such a way that £2 is not adjacent to X\ or £ 3 and £4 is not 
adjacent to £3 or x\. To do this simply choose two adjacent points (labeled x\ and £3) and 
consider the set of all points adjacent to either x\ or £3. There can be at most four such 
points (since the (multi-)graph has degree 3), so we can choose these to be £s,£6,£7,£s- 
Then the remaining points £ 2 and £ 4 cannot be adjacent to X\ or £ 3 . We have thus shown 
that for any topology there exists at least one permutation which will be non- vanishing in 
the light-like limit. Thus the four-point amplitude limit is indeed enough to completely 
determine the correlation function. ^ 

The fact that the four-loop correlation function is compatible at all with the correlation 
function/amplitude duality is itself a highly non-trivial consistency requirement. It would 



14 There is in fact a nice graph theoretical description of whether contributions to completely vanish 
in the light-like limit or not. In order for a contribution from to be non- vanishing in the light-like limit, 
it needs to contain the term 1 / (x^x^x^x^) , cancelling corresponding terms in the prefactor of (|2.23[) . 
If these are not present, the prefactor will force this contribution to vanish. Now any topology will 
contain such a term as long as the corresponding graph (describing f^> and not the numerator polynomial 
contains a J^-cycle (i.e. four points with vertices connected in a square configuration). This is the 
case for I = 4 and £ = 5. 
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be interesting to test this result against the OPE prediction as was done at three loops in 
Section HI 

We can easily continue and do the same analysis at five loops based on the known 
five-loop four-gluon amplitude integrands |4T)] . 

6 Conclusions 

In this paper we have obtained a general expression for the integrand of the three-loop, 
four-point correlation function of stress-tensor multiplets for any gauge group, with four 
undetermined coefficients. In the planar SU(N C ) theory we determine these coefficients 
straightforwardly and uniquely by comparing with the corresponding four-gluon scattering 
amplitude and using the amplitude/correlation function duality. To test this result we 
have compared with the expectations of the OPE in a certain short- distance limit and as a 
byproduct obtained the three-point function (OOJC) at three- loops. We have also obtained 
the four-point correlation function at four loops, albeit without the corresponding OPE 
test. One of the keys to the derivation was the discovery of an unexpected symmetry 
mixing internal and external points in the integrand. 

There are a number of further comments we wish to make. Firstly, in this paper we have 
been mostly interested in deriving correlation functions. In particular, we have emphasised 
the presence of this additional symmetry in the correlation function, where it manifests 
itself as a simple permutation symmetry. However, the techniques and observations in this 
paper have a great deal to say about higher-loop amplitudes as well. Firstly, note that 
a part of the S^+e symmetry remains intact after having taken the light-like limit, and 
thus gives a new hidden symmetry for the four-point amplitude integrand via the ampli- 
tude/correlation function duality. Indeed, this hidden, broken S^+e symmetry explains why 
many of the amplitude graphs appearing in higher-loop four-point amplitudes have equal 
coefficients. It does not, however, fix the overall coefficients. Still, we can get more from 
the amplitude/correlation function duality by observing that the correlation function at £ 
loops contains contributions from products of lower-loop amplitudes as well as the I— loop 
amplitude itself (see for example (13.10P ). Now, if we rewrite the duality relation ( 13. 4 p 
as lim(G4/G4°' ) ) 1/ ' 2 = A±/A^\ we realise that such product terms of non-planar topology 
should not appear on the right-hand side. This gives a consistency condition which fixes 
the coefficients of certain topologies in the correlation function, which in turn then feed 
back in and completely determine many of the coefficients of the £— loop amplitude itself. 
For example, at three loops, the coefficient in front of topology of Fig. 3(b) and hence the 
entire correlation function is completely fixed by (13 . 1 0[) together with the one- and two-loop 
amplitudes and the planarity selection rule we discuss below which eliminates the other 
topologies. Thus the three-loop amplitude can be completely derived from lower-loop am- 
plitudes via the amplitude/correlator duality. At four loops, the coefficients of /a and 
fs in Fig. [10] are completely determined by the contributions of products of lower-loop 
amplitudes, these in turn then determine much (but not all) of the four-loop amplitude. 
Specifically, the coefficient in front of fc is not fixed by this method. On the amplitude 
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side fc contributes to the two amplitude integrals f 2 and d 2 in [35] • Thus, we find we can 
derive all terms in the four-loop amplitude from lower-loop amplitudes, except the (equal) 
coefficients of f 2 and d 2 . Interestingly, the coefficients of all four-loop amplitude integrals 
equal (+1) (these are the coefficients which are fixed by lower loops by this method) ex- 
cept for the f 2 and d 2 coefficients which are then equal to (—1). We currently have no 
explanation for this (see [5T] for an alternative explanation). 

The second comment concerns planarity. We can make an intriguing observation con- 
cerning the topology of the graphs contributing to f^> in the large N c limit. These are the 
graphs given at two loops in in Fig. Upright panel), at three loops in Fig. [3](b) and at four 
loops in Fig. [TOj All these graphs (with the notable exception of one loop) are planar, i.e. 
they can all be drawn on a sheet of paper without any lines crossing, and furthermore all 
other possible contributions to / are from non-planar graphs. If this observation continues 
to higher loops, it would provide a potentially simple and powerful way of greatly reducing 
the number of graphs one need consider in the large N c limit. However it is not immedi- 
ately clear why these graphs for /w should be planar (recall that they do not correspond 
directly to the correlator or the amplitude, but require the attaching of various external 
propagators) . 

What we can argue however is that the graphs contributing to the Born-level (4 + 
tj— point correlation function Gf+ tl (which defines for us the £-loop integrand) should be 
planar in the standard sense. This planar correlation function is related to f^' via (I2.19p . 
Picking a particular component correlation function (for example, choosing the SU(4:) 
channel yf 3 y 24 ) and using (12.171 EES]) ), we find that 

should produce only planar diagrams. Indeed, one finds that insisting on this property at 
three loops uniquely selects topology 3(b) from the four topologies in Fig. [3J The coefficient 
of topology 3(b) will not be determined by this method, but the OPE results in Section 2] 
are enough to fix this remaining coefficient. 

We may thus view the above arguments as an alternative derivation of the three-loop 
correlation function to the one presented in the main text, with the additional benefit that 
this derivation requires no input from the amplitude/correlation function conjecture at all. 
This method however appears to be less useful at higher loops where one is left with more 
and more unfixed coefficients. 

Another comment we would like to make concerns the actual functional dependence of 
the correlation function. What we have done so far is to present the set of Feynman integrals 
which, once computed, will determine the single function of two conformal cross-ratios 
entering in in (I2.23p . The explicit functions for the ladder diagrams in (I4.19p (and for 
the tennis court T in (I4.20p ) are known [48] and they are expressed in terms of standard 
polylogs. It would be very interesting to try to evaluate the two new three-loop integrals E 
and H, Eq. (I4.20p . that we predict. An intriguing question is if some new types of special 
functions will appear there and whether the so-called "maximal transcendentality" rule 
will be respected. We may compare this problem to that of the evaluation of the so-called 
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six-point "remainder function" [ToTfIS] . Recently the result of the Feynman integral two- 
loop calculations has been put in a remarkably simple form [20] involving only classical 
polylogs. We recall that the six-point remainder function depends on three variables, while 
the four-point correlation function depends on two variables. Therefore, it seems that with 
the correlation function we have a better chance to probe this interesting mathematical 
problem at higher loops. 
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A Conventions for the harmonic variables 

Throughout the paper we use two sets of harmonic variables, one for projecting S77(4) 
indices, g\ = (S b ,y b ,) (with the S77(4) indices split as A = (a, a')) and the other for 50(6) 
indices, Y 1 . They are related to each other as follows 

Yi^^i) AB e ab g a A g b B . (A.1) 

Here £/ are building blocks in the expression for the Dirac matrices in six- dimensional 
Euclidean space, see Eq. (B.17) in Ref. [52]. They satisfy 



J2^i) AB ^i) CD = \t ABCD • (A.2) 



7=1 



Using the definition (1 A. 1 [) we find 

(Y 1 ■ Y 2 ) = "£( Y MY2)l = \e ab e a ' b '(y 12 rAy 12 ) b b , = yf 2 , (A.3) 

7 

where (y 12 ) a a , = (yi - y 2 ) a a ,. 

The same matrices {T,j) AB relate the two sets of scalar fields: the six pseudo-real fields 
(j) AB = -(j) BA (with A, B = 1, . . . , 4) and the six real fields $ / (with / = 1, . . . , 6) 

<P AB = ^i) AB ^ , ^ab = \e A BCD<P CD (A.4) 

or equivalently 

Y^ 1 = e ab g\g b B ^ AB . (A.5) 
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In this paper we consider the correlation functions of half-BPS scalar operators 

O ij = tr ($ 7 $ J ) - U IJ ti (® K $ K ) , 

AB,CD = tr ^A B(j) CD^ _ l _ e ABCD tl tyEFfaJ (A 6) 

They are related to each other as 

AB,CD = Ip^ABp^CDQU (A J} 

B Maximal chiral half-BPS nilpotent invariants 

In this appendix we construct the nilpotent superconformal invariants X n appearing in 
the expression for the correlation function (12.191) and explain their properties announced 
in Section I2.2L Similar invariants have been extensively used in the past, also in the 
framework of the Lagrangian insertion procedure, but for four-point correlation functions 
of M = 2 hypermultiplets [3l[53]. There the outer points of the correlation function are 
Grassmann analytic (half-BPS) while the insertion points are chiral (as the M = 2 SYM 
Lagrangian itself). The case M = 4 is rather different because both the outer and insertion 
points are Grassmann analytic. 

A preliminary study of such M = 4 nilpotent invariants was undertaken in [6]. Later on 
a systematic method for building the PSU (2, 2|4) invariants was developed, the supergroup 
formalism of [7HS]- This technique is universal, but it is somewhat technically involved. 
The type of invariant we need here is very special - it has the maximal allowed Grassmann 
degree for a given number of points (see Eq. (I2.17P ) and we are able to present a simple 
alternative construction which makes all properties of the invariant manifest. 

B.l General properties of the M = 4 half-BPS superconformal 
invariants 

The correlation functions (T(l) . . . Tin)) are defined in analytic superspace, which involves 
half of the Grassmann variables, four chiral p a a (see (12. 2p ) and four anti-chiral p",. Part 
of the superconformal group SX(4|4) (the complexification of PSU(2, 2|4)), the chiral 
Poincare supersymmetry Q and the anti-chiral special conformal supersymmetry S act on 
the chiral odd variables p as shifts: 

X ft i jx' ' ft 

d QPa = e a + € a V a' , 

SsP a a =<(e a + liy a a), (B.l) 

and analogously for the variables p and generators Q and S. The total number of odd 
parameters e,^ in (IB . 1 j) is 16. Using this freedom, we can gauge away four sets of analytic 
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odd variables, for example, 



Q + S gauge: p x = p 2 = p 3 = Pa = , (B.2) 

as well aJ^I 

Q + S gauge: p x = p 2 = p 3 = h = . (B.3) 

If we restrict ourselves to the four-point correlation function (T(l) . . . T(4)), then the entire 
dependence on the odd variables is completely fixed by superconformal symmetry: given 
the p = p = component (12. 9p . the global supersymmetry transformation which leads to 
the gauges (IB.2I) . (IB.3I) unambiguously restores the dependence on p, p. 

However, with n > 5 points we can have nilpotent superconformal invariants. Indeed, 
the gauges (IB.2I) . (1B.3|) completely exhaust the power of the odd part of SX(4|4); any 
remaining p^ (or pi) with % = 5,...,n will automatically be supersymmetric invariant. 
Further, the Z4 centre of the R symmetry group SU(4) acts on the odd variables by 
rescaling: p — >■ up and p — > u~ l p with u 4 = 1. Thus, we have the following types of Z4 
invariants: chiral (p) 4fc , anti-chiral (p) Ak and mixed {pp) h , or combinations thereof. In each 
case the Grassmann degree cannot exceed the value An — 16 corresponding to the number 
of odd variables left at our disposal after fixing the gauges (1B.2P and (IB.3p . We shall 
call maximal the nilpotent invariants which involve all the available odd variables. For 
instance, the maximal chiral invariant is of the form {p) A ^ n ~ A \ Notice that such invariants 
are unique, up to an overall factor depending only on x and y. The dependence on the 
harmonic variables y is further restricted by harmonic analyticity (i.e., irreducibility under 
SU(4)), which allows only polynomial dependence of a fixed degree, in our case of degree 
four (matching the U(l) charge of the stress-tensor multiplet T). 

Now, the subject of this paper are the four-point correlation functions (T(l) . . . T(4)). 
As explained above, their dependence on the odd variables is obtained by transforming the 
bosonic variables x, y under conformal supersymmetry with the parameters which have 
been used to achieve the gauges (1B.2|) . (1B.3H . Such transformations always create pairs 
Pipi. Therefore, if we decide to eliminate the anti-chiral variables by fixing the gauge (IB. 3D . 
we automatically eliminate the chiral ones as well. 

The situation is different for n > 5: There, even in the gauges (IB.2p . (IB. 31) . we can have 



purely chiral (and anti-chiral) or mixed nilpotent invariants. But again in the context of 
this paper, for the purpose of Lagrangian insertions into the four-point correlation function 
(see (12.141) ). the relevant part of the higher-point correlation functions are the components 
in which the M = 4 SYM chiral Lagrangian appears at n — 4 points. This Lagrangian 
belongs to the purely chiral sector of the stress-tensor multiplet, see (12. 6p . For this reason 
we may ignore the p dependence of the correlation function, both for n = 4 (when it is 
fixed by supersymmetry) and n > 4 (when it is irrelevant for our purposes). This justifies 
our decision to set all pi = 0, with % — 1, . . . , n. 



15 In the complexified description p and p are treated as independent, so the two gauges can be fixed 
separately. 
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So, according to the insertion formula ( I2.14p . we are interested in the n— point corre- 
lation function (I2.13p . Gf^ e . e ~ (ps) 4 • • • (Pa+i) 4 + ■■■■, where the dots denote the super- 
conformal completion to a full n— point purely chiral nilpotent invariant of the maximal 
type. Consequently, its dependence on the odd variables is uniquely determined by super- 
symmetry. In the next subsection we construct such an invariant and show that the only 
freedom in it is a factor depending only on x. Most importantly, our construction makes 
the permutation symmetry S n with respect to all points manifest. 

B.2 Constructing the maximal nilpotent invariant 

Consider the following chiral Grassmann integral 

/n 
d\ d'e'd^d^ 1 J] 6^ (p& - (f a + ety*) - < + gy&) ) . 
i=l 

(B.4) 

It depends on n points (xi,yi,pi) in the chiral analytic superspace and has the maximal 
Grassmann degree 4n — 16. The invariance of the integral under shifts of the integration 
variables clearly implies invariance of X n under Q and S supersymmetry (IB. II) . 

Let us examine the properties of X n under the conformal inversion X{ — > x" 1 , — > x~ x po 

Z n {x-\y l ,x- 1 p l )= I d 4 ed 4 e'd 4 td^' \{{xj)- 2 ^ ( Pi - x t {e + e' Vi ) - (I + ? Vi )) . 

i 

Exchanging the integration variables e -H- ^ and e' -H- we arrive at 

n 

In^iViiX^Pi) = Yl ( X k)~ X 1n{.Xi,yi,Pi) ■ (B.5) 
k=l 

In addition, it is easy to show that X n is also invariant under shifts Xi — > Xi+5. Since special 
conformal transformations are a combination of translations and inversions, K = IPI, we 
see that X n (xi,yi, P i) is conformally covariant with weights (—2) at each point. 

The R symmetry group SU(A) can be treated in complete analogy with the confor- 
mal group SU(2, 2) (both being real forms of SL(4, C)). Its action can be reduced to a 
combination of translations of y a a , and inversion in the group space, 

y a a > -> {y'X , P a ->P a (rt'- (B.6) 

We verify that X n is covariant under SU (4) as well, 

n 

ln(xi, y' 1 , pi y' 1 ) = Yl {ylY x In(xi, yi, pi) . (B.7) 

k=l 

We remark that, being a chiral object, X n (xi,y.i, pi) is not invariant under the anti- 
chiral part of Poincare supersymmetry which acts only on the x's, 5qX aa = p a e a . The 
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same applies to the chiral half of special conformal supersymmetry S = IQI. In fact, in 
the gauge (1B.3[) all the parameters e and £ are fixed and Q + S supersymmetry is completely 
broken. If needed, the dependence on p and with it the Q + S invariance may be restored 
by undoing the gauge fixing, but this is not necessary here. 

To summarise, the nilpotent object X n constructed in (1B.4I) is: 



• chiral and Grassmann analytic of maximal Grassmann degree 4(n — 4); 

• invariant under Q + S supersymmetry; 

• covariant under SU(2,2) x £77(4) with conformal weight (—2) and SU(A) weight 
(+4) at each point; 

• harmonic analytic, i.e. polynomial in the y's of degree four in each matching the 
£[7 (4) weight; 

• fully symmetric under the exchange of any pair of superspace points (xi,yi, pi). 

What is the freedom left in this invariant? As explained earlier, the p— dependence of 
X n is already fixed by supersymmetry. We could modify X n by multiplying it by a bosonic 
function f(xi,yi). For our purposes we need an invariant with £77(4) weights (+4) at each 
point, to match the properties of the stress-tensor multiplet T. The invariant X n already 
has the required weight, therefore the function f(xi,yi) must depend on harmonic cross- 
ratios (the y— analogs of the space-time cross-ratios (I4.18P ). However, the latter inevitably 
introduce harmonic singularities (poles in yfj) which are forbidden by harmonic analyticity. 
We conclude that the only freedom left in the invariant X n is a function f{xj). It must 
carry conformal weight (+4) at each point in order for X n (xi, Pi, yi) x f(xi) to match the 
weight (+2) of the T's. Finally, the complete crossing symmetry of the correlation function 
(T(l) . . .T(n)), combined with the point permutation symmetry of X n , imply that f(xi) 
must have full symmetry £„. This proves our claim made in f)2.19p . 

B.3 Gauge- fixed form 

Evaluating the Grassmann integrals in (1B.4[) in general is cumbersome and leads to a 
complicated Jacobian involving x and y. The calculation is greatly facilitated by choosing 
the (Q + £) gauge (IB. 21) . p\ — . . . — p± — 0. In addition, we may use part of the conformal 
and £77(4) symmetry (translations and inversion) to fix the 

££7(2,2) x £7/(4) gauge: x 2 , y 2 ^ oo , x 3 = y 3 = . (B.8) 

The rest of SU (2, 2) x ££7(4) allows us to bring the remaining x and y variables to the 
diagonal 2x2 matrix form P|[TU] 



XuX t 1 



z 
z 



2/i42/i 1 



w 
w 



(B.9) 
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In the gauge (lB.8p . flB.9j) the conformal cross-ratios (I4.18P becom 



u 



xA 



(l-z){l-z) 



X 



and similarly for the y— variables: 



2 2 — = (1 - W)(l - , 



2 ^ > 



2/232/41 . 2/4I 

-5 — 5- — > = WW . 



(B.10) 



(B.n; 



Then, the product of the first four delta functions on the right-hand side of flB.4|) reduces 
in the gauge flB.8|) to 



(xlylxtyirS^\^\e)5^\^(e'y 41 y^ - x^e') 



(B.12) 



so that the integration over e, £, £' in flB.4|) is trivial. To integrate over e' we introduce a 
parametrisation for its components and use ( IB.lOj) and (IB. lip to recast the argument of 
the last delta function in (1B.12j) in the form 



a (3 
7 5 



a(w — z) (5{w — z) 
j(w — z) S(w — z) 



(B.13) 



Thus, in the gauge ( IB.lOj) and ( IB. lip the nilpotent invariant (1B.4|) becomes 



X„| Pl= ... =P4= o -> (x\xlylyl) (w - z)(w - z)(w - z){w - z) x Jj5 (4) (pi) . (B.14) 



i=5 



Let us now examine the expression for the rational factor R from (12.181) in the same gauge. 
It is easy to check using (IB. 10p and (IB. 1 1 [) that it takes the following form [9] 



(xlsxl^xl^x^) x #(1,2,3,4) -»■ (x\xlylyl) 2 {w - z)(w - z)(w -«)(«;- z) 



Comparison of (IB.14j) and (IB. 15[) yields 

n 

3-n\pi=...=p4=0 = { x lZ X 2i X 12 X 2'i X ?,A X Al) X -^(1) 2, 3, 4) X J |^ <5 ^ \pi) , 

exactly as announced in (j2.17|) . 



(B.15) 



Ti.16) 



t=5 



16 The variables z, z are real and independent from each other in the case of Minkowski metric, and 
complex conjugate to each other in the Euclidean case [10] . Similar variables have been used in the past 
for solving the superconformal Ward identities for four-point correlation functions {51151-HT] . 



39 



C Numerical evaluation of the short-distance limit 



In this appendix, we study the asymptotic expansion of the Feynman integrals defined in 
f)3.14p and (I3.13P in the Euclidean region x\ — > 12, 13 ^ ^4- In terms of the conformal 
cross-ratios ( I4.18p . this limit corresponds to u — >■ 0, v — > 1. As was explained in Section 4, 
the advantage of this limit compared to the conventional short-distance limit X\ — >■ x 2 , 
corresponding to u — > with v fixed, is that it allows us to single out the contribution to 
the OPE (14.11) of a couple of operators with the lowest spin and scaling dimension. 

We start with the one-loop ladder function g(l,2,3,4) and examine the 'three-star' 
integral in (14.281) . At the first step, we dimensionally regularise the integral and rewrite it 
using the standard Feynman parameter technique: 

d4 ~ 2€x5 ^(2 - e) r(i + e) f , dqid f q3 6(1 - qi ~ q2 ~ q t • (ci) 



4vr 2 J x 2 15 xj 5 x 2 35 4 J (qiq 2 x 2 12 + qiq 3 x 2 13 + q 2 q 3 xl 3 ) 1+e 

Then, we follow the usual Mellin-Barnes (MB) procedure for rewriting the sum of terms in 
the denominator as a product at the cost of introducing two new (contour) integrals with 
variables z\ and z 2 (with fixed real part and varying imaginary part) 



x [\q 1 dq 2 dq 3 5(l-q 1 -q 2 -q 3 )q; {1+e+Zl \~ {1+e+Z2) q^ +z \ (C.2) 
J 

where the integration contours run along the imaginary axis to the left of the poles of the 
first two T— function and to the right of the poles of the third T— function. Finally we 
perform the Feynman parameter integration (sec for example [49J for more details) to get 

r(2-e) [ +lca d Zl dz 2 



4r(l - 2e)(40 



J -ioo 



.22 / 2 \ z l 
X Y1 \ I X 23 



xT(l + z 1 + z 2 )T(l + e + z 1 + z 2 ) -fM . (C.3) 

\ x 13/ V x 13/ 

At first sight, it may seem surprising that we use a dimensionally regularised form of the 
integral (ICIO . given the fact that it is finite for e — > and hence conformally covariant. The 
reason for this can be understood in more complicated integrals, for example, the three- 
loop ladder (see, e.g., [IE]). The MB parametrisation of the integral produces factors of 
r(— 2e) in the denominator, which are compensated by poles arising due to the integration 
over the Zi variables. The result is of course finite and non-vanishing. The fact that this 
standard technique works in D = 4 — 2e does not clash with conformal invariance precisely 
because the integrals are finite: We reduce a conformal integral in four dimensions to a 
non-conformal three-point integral. Although the latter is computed in D = 4 — 2e, in the 
resulting finite expression we can simply upgrade the 0(e°) part of the representation to a 
conformal four-point integral. The procedure leads to a substantially reduced number of 
MB parameters as opposed to what one would obtain starting from the four-point integrals. 
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Coming back to our example (ICIO . the integral is finite as anticipated and we can 
safely put e = 0, leading together with (14.281) to 

g(l, 2, 3, 4) = / y^f ^(-^(-^(1 + Zl + z 2 )f «« . (C.4) 

The MB.m package [51] may be used to evaluate this integral. Here it suggests we put 
the real parts of the integration contours to {—3/16, —3/8} for {2:1,2:2}, respectively. For 
u — > and v — > 1, we obtained with the help of the MBasymptotics.m package [55] 

0(1,2,3,4) = Qln«-i + 0(«)V (C.5) 

To match our Wick rotation conventions (positive Euclidean signature) the output of the 
Mellin-Barnes software has to be changed by a minus sign per loop order. 

For the two-loop integral h, Eq. ( I3.13p . and the three- loop integrals L, E, H, Eq. (13.141) . 
there is no complete point-exchange symmetry, so each of the integrals entering (I4.26P must 
be considered separately. To obtain an MB representation for h(l, 2; 3, 4) we send point 4 
to infinity upon which either subintegral becomes a three-star like (IC.lj) . yielding two MB 
parameters each. For the three-loop ladder L and the easy three-loop integral E we start 
with the outer subintegrals, the middle integral is then a three-star as before so that we 
obtain representations with a total of six parameters (see (13.141) for the integral definitions). 

The hard integral H is much more nested. Let us send point 4 to infinity as before. 
We would prefer to start with the central integral which becomes a three star, aiming at a 
9-parameter representation. Unfortunately, the MB.m package fails to find an integration 
contour for this representation!^! Starting with the outer integrals and leaving the central 
integration to the end produces a 12-parameter representation which we can run; its 0(e°) 
part contains several integrals with up to nine parameters. 

In the asymptotics the number of parameters is smallest in the leading log part; as 
observed also on pole parts of IR or UV divergent integrals the dimensionality of the 
MB integrals increases towards the finite part. Yet we can obtain rather good numerical 
precision because we compute only number coefficients and not functions. After (for the 
hard integral rather heavy) application of the barnesroutines.m package [55], the worst 
parameter count occurs for the finite part of the hard integral, which in the u — > 0, v — > 1 
limit contains a few five-parameter integrals. For up to four parameters the Cuhre routines 
in the Cuba library [56] allow excellent precision, but in this case we could not obtain more 
than three significant digits. 

The coefficients that we find are either simple rationals like 3/2, 5/4 or simple rationals 
times C(3) or C(5), and linear combinations thereof. The denominators of the coefficients 
of the subleading In u divergences are powers of 2 wherever the MBasymptotics.m package 
leads to analytic results. Using this as a guideline it is very easy to guess exact results from 

17 The same difficulty was encountered when studying the pentabox in [31] . The hard and easy inte- 
grals both have a pentabox subintegral; we were in fact surprised that the problem did not occur when 
constructing the representation for the easy integral E described above. 



41 



the usually excellently precise numerical output as long as there is no mixture of rational 
numbers and values. 

In this way, we find the following results for asymptotic behaviour of loop integrals 
entering (14.261) as u — > and v — > 1 (we display the leading terms onlyg 
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— (Inn) 2 - 


3 , 

— mu + 
1 p, 




x^L(l,2;3,4) = 
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16 V ; 



2 5 . 5 

H In it 

32 16 

2 5 , 
H In u 

32 

5 

32 ' 



9 17 
— lnn+-C(3)- — 
32 4 SV ; 16 

(C.6) 

Using these formulae we can find the short-distance limit of the loop corrections (12.221) to 
the four-point correlation function displayed in (I4.3ip . 
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